Nonanalytic paramagnetic response of itinerant fermions away and near a 
ferromagnetic quantum phase transition 
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We study nonanalytic paramagnetic response of an interacting Fermi system both away and in 
the vicinity of a ferromagnetic quantum phase transition (QCP). Previous studies found that (i) 
the spin susceptibility x scales linearly with either the temperature T or magnetic field H in the 
weak-coupling regime; (ii) the interaction in the Cooper channel affects this scaling via logarithmic 
renormalization of prefactors of the T, \H\ terms, and may even reverse the signs of these terms 
at low enough energies. We show that Cooper renormalization becomes effective only at very low 
energies, which get even smaller near a QCP. However, even in the absence of such renormalization, 
generic (non-Cooper) higher-order processes may also inverse the sign of T, \H\ scaling. We derive 
the thermodynamic potential as a function of magnetization and show that it contains, in addition 
to regular terms, a non-analytic \M\^ term, which becomes M"^ /T at finite T. We show that 
regular (Af'^,M*, ...) terms originate from fermions with energies of order of the bandwidth, while 
the non-analytic term comes from low-energy fermions. We consider the vicinity of a ferromagnetic 
QCP by generalizing the Eliashberg treatment of the spin-fermion model to finite magnetic field, 
and show that the \M\^ term crosses over to a non- Fermi- liquid form \M\''^^ near a QCP. The 
prefactor of the \M\'^'''^ term is negative, which indicates that the system undergoes a first-order 
rather than a continuous transition to ferromagnetism. We compare two scenarios of the breakdown 
of a continuous QCP: a first-order instability and a spiral phase; the latter may arise from the 
nonanalytic dependence of x on the momentum. In a model with a long-range interaction in the 
spin channel, we show that the first-order transition occurs before the spiral instability. 

PACS numbers: 71.10. Ay, 71.10 Pm 



I. INTRODUCTION 



The Landau Fermi-liquid (FL) theory postulates that 
at low enough energies a system of interacting fermions 
behaves as a weakly-interacting gas of quasiparticles with 
renormalized parameters: effective mass, Lande g-factor, 
etCfi The thermodynamics of a canonical FL is con- 
structed under the assumption that the fermion-fermion 
interaction is absorbed entirely into a set of the renormal- 
ization factors (Landau parameters), while the residual 
interaction between quasiparticles can be neglected. In 
this approximation, the FL behaves as a Fermi gas of 
free quasiparticles. In particular, the specific heat coef- 
ficient, j{T,H) = C{T,H)/T, and the uniform, static 
spin susceptibility, x(T, 7?), remain finite in the limit of 
T,H — > 0, while their T and H dependences follow the 
familiar Sommerfeld expansions in powers of and H^. 

It has long been known that neglecting the residual 
interaction leaves some important physics behind. In 
particular, non-trivial kinetics of a FL is entirely due 
to the residual interaction among quasiparticles. The 
effect of the residual interaction on thermodynamics of 
Fermi systems has been studied intensively in recent 
years (for a review, see Refs. 33). It is well estab- 
lished by now that T and H dependences of j{T, H) 
and x(r, i?) are nonanalytic. In two dimensions (2D), 
both 7 and v are linear rather then quadratic in T and 
\H\ (Refs. iiil31iill31illlllllllii0l3li30). in ad- 
dition, the nonuniform spin susceptibility, x{q)i scales 



linearly with \q\ for q-^kp (Refs. ItHTqI). 

The nonanalytic behavior originates from a dynamic, 
long-range component of the residual interaction medi- 
ated by virtual particle- hole pairs. Two regions in the 
space of momentum transfers contribute to the long- 
range dynamics. The first one is the region of small q, 
where the long-range interaction arises due to the fJ/g 
form of the fermion polarizability (this form is also the 
reason for Landau damping). In real space, this com- 
ponent of the interaction falls off slowly, e.g., as Vl/r 
in 2D. The second one is the region around 2kF, where 
the Kohn anomaly generates dynamic Friedel oscillations 
falling off as ncos{2kFr)/r^/'^ in 2D. With this in mind, 
the non-analytic behavior of the free energy can be ob- 
tained by the following scaling argument. The range of 
the interaction via particle-hole pairs is determined by 
a characteristic size of the pair, Lph, which is large at 
small energy scales. At finite temperature, Lph ~ vp/T 
by the uncertainty principle. To second order in the 
bare interaction, two quasiparticles interact via a sin- 
gle particle-hole pair. The energy of order T, carried 
by such a pair, is distributed over a volume L^j^ oc T^^ . 
The contribution from such process to the free energy per 
unit volume is of order 5F ~ v?T / L^^i °c T^^"^ , where 
u is the dimensionless coupling constant. Consequently, 
7(T) = -d'^F/dT^ cx T^"^. Likewise, at T = but in 
finite magnetic field, a characteristic energy scale is the 
Zeeman splitting 2fj,B\H\ and Lph ~ vp / fJ-B\H\. Hence, 
6F cx \H\^^^ and x{H) (x \H\^^\ For D ^ 2, this 
implies that 7(T) cx T and x{H) cx \H\. For D = 3, 
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power-counting misses logarithmic factors which are re- 
covered by an explicit calculation. 

A perturbation theory indeed shows that 7 and x 
depend linearly on T and in 2D, and as InT and 
hi\H\ in 3D. To second order in the interaction, it has 
been iounS'^^''^^ that 

57(T, H = Q) = [/f (tt) + 3/f (vr)] ^7o'fla) 

5x(r = 0,i/) = /,^(7r)^Xr (lb) 
6x{T,H^Q) = fA^)^xr (Ic) 

Sx{T^O,H = 0,q)^^f^{n)^xr, (Id) 

where (tt) = (to/tt) [[/ (0) - C/ (2fcF) /2] and /.(vr) = 
— {m/2TT)U{2kp) are the charge/spin component of the 
(first-order) backscattering amplitude fc/s{(^ — & 
is the angle between the incoming momenta of two 
fermions. Also in Eqs. (|lamd|) . 7q^ = rmr/3 is the spe- 
cific heat coefficient of a 2D Fermi gas, Xo^ = A's'tt-/'"' is 
the spin susceptibility of a 2D Fermi gas, e_F is the Fermi 
energy, is the Bohr magneton, and all relevant en- 
ergy scales-T, ^sji?!, and vp |q|-are small compared to 
ep- (The scaling forms as functions of all three variables 
can also be obtained, see Ref. [ll|.) Scattering processes 
contributing to Eqs. (Ilajlldp are characterized by spe- 
cial kinematics ("backscattering"): two fermions move 
in almost opposite directions before a collision and then 
either continue to move along the same path (momen- 
tum transfer q = 0) or scatter back (momentum transfer 
2kp). 

The intriguing feature of the perturbative results is 
that the spin susceptibility is not only nonanalytic but 
also an increasing function of all three arguments: 
T, and q. Since one should expect the susceptibility to 
decrease at least at energies much larger than ep, a. nat- 
ural conclusion is that x a maximum at interme- 
diate energies. If this behavior survives beyond weak- 
coupling, it implies non-trivial consequences for a mag- 
netic phase transition in such a system. Indeed, a max- 
imum of x(2^7 H,q — 0) at finite H gives rise to a lo- 
cal minimum in the free energy at finite magnetization 
M. As x{^ = 0) increases, this minimum becomes de- 
generate in energy with a non-magnetic state implying 
that a ferromagnetic state emerges via a discontinuous, 
first-order transition accompanied by a metamagnetic re- 
sponse away from the critical point. On the other hand, 
a maximum of x{T, H = 0,q) at finite q implies that 
the system may also undergo a transition into a spi- 
ral rather than uniform magnetic state. Both scenarios 
imply a breakdown of the Hertz-Millis-Moriya (HMM) 
model of a continuous, quantum, ferromagnetic phase 
transition22ii^i2^. The first-order instability has been dis- 
cussed in recent literature.— 1^ It is not clear, however, 
which of the two instabilities-the first-order or spiral one- 
occurs first. One of the aims of this paper is to clarify 
this issue. 



Experimentally, a linear T dependence of the specific 
heat coefficient has been observed in monolayers of '^He 
(Ref. [[2^); both the sign and the magnitude of the ef- 
fect are consistent with Eq. (ITa| (Ref. s''^. For the spin 
susceptibility, the experimental situation is less clear. A 
quasi-linear dependence of x on T was observed in a Si- 
based 2D heterostructure^'*; however, the slope is oppo- 
site in sign to that in (jlcp . On the other hand, a num- 
ber of experiments on this and other heterostructures (Si 
Ref. m, n-GaAs Ref. M, and AlAs Ref. [13) have found 
that X increases with magnetization, in agreement with 
Eq. (fTb]) . A linear temperature dependence of x has 
recently been observed in the normal phase of Fe-based 
pnictidesj^ the sign of the slope is consistent with Eq. 

A finear Igj-dependence of x{T — 0,H — 0, g) has 
recently been proposed to infiuence ordering of nuclear 
spins via a Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction mediated by interacting rather than free 
electrons Because of the \q\ term, the dispersion of 
nuclear spin waves in the RKKY-ordered state, ujs{q) oc 
xil), is linear rather than quadratic in q. In 2D, this 
implies that the nuclear magnetic order is stable with re- 
spect to thermal fluctuations, which opens a possibility 
to freeze nuclear spins at experimentally accessible tem- 
peratures with potential applications in quantum com- 
puting. 

Conflicting observations of the temperature and 
magnetic- field dependences of xiT,H) and potential ap- 
plications in quantum computing call for a detailed the- 
ory of the nonanalytic effects in the spin response of 2D 
and 3D Fermi systems. In particular, it is important to 
understand whether the weak-coupling results can be ex- 
tended into a non-perturbative regime near a ferromag- 
netic transition. 

Several groups have recently investigated this 
issueJ^iHiiiii^ii^iiiii^i^i^ It turns out that the result 
for the specific heat is robust: for D < 3, all higher- 
order corrections can be absorbed into renormalization 
of the backscattering amplitudes fc (tt) and (tt) i n 
the second-order resuh, Eq. ^ (Ref. [Ullllilll) . 
One particular consequence of this result, which still 
awaits for an experimental verification, is the addi- 
tional logarithmic dependence of the specific heat 
coefficient 7(1") oc T^^^ /(Inep/T)^ resulting from 
renormalizations of /^/s (tt) in the Cooper chan- 
nel (4^(7r) cx \n-^{ep/T) in the limit of T 0, 
Refs. [ijllllll). In 3D, there are additional T^lnT 
nerms in 7, which are not expressed via backscattering 

13.33 

For the spin susceptibility, the situation is more com- 
plex: even in 2D, not all higher-order processes can 
be absorbed into renormalization of the backscattering 
amplitudes in the second-order results. The remaining 
processes do not have special kinematics: the momenta 
of incoming fermions are not correlated and momenta 
transfers are generic rather than peaked either near or 
near 2kp. The signs of these extra linear contributions 



3 



to x(2^: -f^) alternate with order of the perturbation the- 
ory, which opens a possibiUty for sign of x(2^i H) to be 
reversed upon resummation. In addition, the backscat- 
tering contribution is suppressed by Cooper logarithms, 
leaving the non-backscattering processes as the main con- 
tributors to linear in T and \H\ terms in the susceptibility 
at sufficiently low 

In this paper, we develop a general theory of the non- 
analytic behavior of the spin susceptibility in two and 
three dimensions, both in the FL-regime and also in the 
vicinity of a ferromagnetic quantum critical point (QCP). 
In Sec. ini we discuss the 2D case. After a brief re- 
view of the perturbation theory for x in Sec. Ill Al we 
construct in Sec. Ill Bl an expansion of the exact suscep- 
tibility in skeleton diagrams with an increasing number 
of dynamic polarization bubbles. Physically, such an ex- 
pansion corresponds to collecting all processes involving a 
given number of virtual particle- hole pairs. In Sec. Ill B II 
we show that all diagrams with two dynamic bubbles give 
effectively second-order results (|lbllcp but with the exact 
rather than perturbative backscattering amplitudes. In 
Sec lII B 2l we consider processes with more than two dy- 
namic bubbles and show that they also give rise to linear 
T and \H\ terms in x{T,H). We evaluate the diagrams 
up to fourth order in dynamic bubbles and calculate x 
explicitly for a model form of the scattering amplitude 
parameterized by the ffist two harmonics, fs^ and fs^i- 
In Sec. Ill CI we address an issue of the sign of the T- and 
i?-dependences of x- We show that higher-order process 
can reverse the sign of backscattering contribution for a 
strong enough interaction, even if logarithmic renormal- 
izations in the Cooper channel are neglected. In the same 
Section, we also analyze the role of Cooper renormaliza- 
tions for a system with a short-range interaction and for 
a 2D electron gas with Coulomb interaction in the large 
iV limit, relevant mostly for valley-degenerate semicon- 
ductor heterostructures. In agreement with Ref. [H, we 
find that the slope of x{T, H) in a Coulomb gas changes 
sign below a certain energy; however, this energy is of 
order E* = ep exp {—N^/'^ /^/2) in the large-A^ model. 
Already for the case of two valleys (TV = 4), E* is too 
low for this mechanism to be responsible for the observed 
negative sign of the T-dependence of x in s Si-MOSFET 
(Ref. 24). 

Next, we obtain a general form of the thermodynamic 
potential for a 2D FL with an arbitrary strong interaction 
(Sec. HID 1[) and extend the analysis of the magnetic-field 
and temperature dependences of x(r, H) to both FL- and 
non-FL regions near a ferromagnetic QCP in 2D (Sec. 
imp . In Sec. nil D 21 we neglect Cooper renormalizations 
and show that while x(2^i H) increases with iJ, T in both 
regimes, the \H\,T scaling holds only up to a certain 
energy which decreases as the QCP is approached. At 
higher energies, the magnetic-field and temperature de- 
pendences of x(T,H) are \H\^/-^ and TlnT, respectively. 
The increase of x with H signals an imminent breakdown 
of the continuous ferromagnetic transition. We discuss 
possible scenarios of quantum- and finite-temperature 



ferromagnetic phase transitions in Sec. IIIIEI In partic- 
ular, we show that for a large radius of the interaction 
in the spin channel the first-order transition always pre- 
empts the spiral instability. Finally, in Sec. IIIIFI we 
show that the increase of x with H and T near a QCP 
is not affected by renormalization in the Cooper chan- 
nel, as this renormalization becomes relevant only below 
an energy which decreases exponentially as the QCP is 
approached. 

In Sec. IIV|, we consider the 3D case. In Sec. II V H 
we show that cx H^\n\H\ form of x{T — 0,H) in a 
3D FL transforms into a weaker, if^lnln|i?| form near 
a ferromagnetic QCP. In Sec. II V 21 we analyze the T- 
dependence of x{T, H — 0) in 3D. A 3D FL is peculiar in 
a sense that xC^) scales as without an extra logarith- 
mic factor— i^^. We generalize the earlier result for the 
scaling by Beal-Monod et al.— and show that the pref- 
actor of the term is non-universal: its magnitude and 
sign depend on details of the fermion dispersion. For the 

dispersion discussed in Ref. HI, the prefactor of the 
term is negative, i.e., xiT) decreases with T. However, 
X may increase with T for a more complex dispersion. A 
increase of xiT) with T has been observed in a number 
of exchange-enhanced paramagnetic metals.^ 

Section |V] summarizes our conclusions. Some technical 
details of the derivations are given in Appendices [A]- [C] 
Some of the results presented here were published in a 
shorter form in Ref. [ij. 



II. MAGNETIC RESPONSE OF A 2D FERMI 
LIQUID 

The spin susceptibility at zero temperature and in zero 
magnetic field, x{T = 0, iJ = 0) is described by the con- 
ventional FL theory (Ref.[ll). The subject of our study is 
the temperature- and field-dependent part of the suscep- 
tibility: Sx{T,H) = x{T,H)-x{0,0). The most straight- 
forward way to obtain 5xiT, H) is to evaluate the thermo- 
dynamic potential S (T, H) and differentiate it twice with 
respect to the field. In contrast to the linear-response 
theory, which generates a large number of diagrams, the 
number of relevant diagrams for the thermodynamic po- 
tential is rather small. The prefactor of the term in 
the thermodynamic potential gives the T-dependent spin 
susceptibility, while the nonanalytic \H\^ term gives the 
field-dependent (nonlinear) susceptibility. 

A. Second-order perturbation theory 

To second order in the interaction 5x{T, H) was con- 
sidered in RefsHidH where it was found that 

• x(^i l-^l) is nonanalytic in both arguments and 
scales as max{T, |i?|}); 

• the nonanalyticity comes from the states near the 
Fermi surface; 
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• only 2fci? scattering is relevant, thus the prefactors 
of the linear terms in T and in H contain only the 
2kp component of the unteraction U{q). 

In this Section, we overview briefly the second-order per- 
turbation theory, because later we will need to under- 
stand what replaces U{2kp) in the interaction vertices 
beyond the second order. 

At second order in U{q), the field-dependent part of 
the thermodynamic potential 'E.{T,H) is given by a sin- 
gle diagram shown in Fig.[TJ In this diagram, the spins of 
fermions in one of the bubbles are opposite to those in an- 
other bubble. The nonanalytic contribution to S(r, 
originates from the 2kp nonanalyticity of the dynamic 
polarization bubble in zero field and, hence, is propor- 
tional to U{q — 2kp). Finite magnetic field cuts off the 
nonanalyticity, but at a price that the derivatives with 
respect to the field become nonanalytic in H. With all 
four fermions near the Fermi surface, the 2kp diagram 
necessarily contains two spin-up fermions with momenta 
near hp and — k^^ and two spin-down fermions also with 
momenta near hp and —hp. These four fermions can 
be re-grouped into two up-down bubbles, each with a 
small momentum transfer. This simplifies the computa- 
tions substantially because the polarization bubble has a 
much simpler form for small q than for q near 2kp. 

The magnetic field enters the problem via the Zeeman 
shifts of single-fermion energies in the Green's functions 



FIG. 1: The field-dependent part of the thermodynamic po- 
tential at second order. k,p,q are the four- momenta: k = 
(k,Lj,„), etc. 




The nonanalytic part of S2 (T, H) is obtained by keeping 
the square of the dynamic term in Eq.®, i.e., replacing 
n^j^ by J^^Pfi- This gives 



S2 {T,H)^ 



■T 



'■2k 



^rVrii In 



^^n 

{^l^-l^Y^{vpqy 



(8) 



where 



1 



iijjm - Ck ± A/2' 



where 



A EE 2^B-ff 



(2) 



(3) 



is the Zeeman energy. The up-down polarization bubble 
is defined as 

(g, a«) = Gt (w™ + f7„„ k + q) (w™, k) , 

k 

(4) 

where T^^, is a shorthand for / d^fc/(27r)^. For 

q <^ kp, the up-down bubble can be separated into the 
static and dynamic parts as 



iin{q,nm)^-iy{i-Pti), 



(5) 



where v = m/2TT is the density of states at the Fermi 
surface and 



(6) 



Re-expressed in terms of the up-down bubbles, the dia- 
gram in Fig. [1] reads 



S2 (T,iJ) = - 



?72(2fc/ 



^r^n?^(f^„„g). (7) 



U2kp = vU{2kp) 



(9) 



and W is the high-energy cutoff which, in general, is of 
order of the bandwidth. 

The logarithm in the frequency sum in Eq. ([8|) orig- 
inates from the f^m/kl form of the polarization bubble 
at vpq 3> rim, A, i-e., from the long-range tail of the 
dynamic bubble (in real space, fim/l^l transforms into 
^Im/f)- If not for the logarithm, 'E.2{T,H) would be ex- 
pandable in powers of and H^. The logarithm breaks 
analyticity. Replacing the Matsubara sum by a contour 
integral, and subtracting off the field-independent part, 
we obtain from Eq. ^ 



S2 (T,i/) 



|A| 



diin-^ coth 



n \ 

2t) 



(10) 



The integral in Eq. (fTU|) can be solved exactly (in terms of 
polylogarithmic functions), but we actually do not need 
this solution, as the T- and iJ-dependent spin suscepti- 
bility can be obtained directly from Eq. (|10p by differen- 
tiating it twice with respect to H. This yields 



Sx2{T,H) = - 



2 2 

= u 



— 51 A 

^''^ 2Ep \ 2T 



(11) 



where ~ fi^m/TT is the spin susceptibility of a free 
2D Fermi gas, and the scaling function S{x) is 



S{x) — coth a; 



2 sinh X 



(12) 
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FIG. 2: Diagrams for the thermodynamic potential beyond 
second order. 




The asymptotic limits of S are S{x — > 0) = l/2x and 
S{x ^ oo) ^ 1. Substituting these limits into Eq. pT|) . 
we find that the susceptibility increases linearly with the 
largest of the two energy scales, T and A. Schematically, 

6x2{T,H) = ul,^^xl'' (13) 

where E = maxjT, |A|}. 

If the same calculation is performed in real rather 
than Matsubara frequencies, the frequency integral con- 
tains the product of the real and imaginary parts of 
the retarded, dynamic bubble: U2i^^KeP^^lmP^^. This 
allows for a transparent physical interpretation of the 
two-bubble diagramji^ Indeed, Imll^ can be thought 
of as the spectral density of particle-hole pairs, while 
U2j,^Ren^ as of the dynamic interaction between the 
two fermions in the particle-hole pair. The product 
dn(fqul^.^ReU"^luiU^' cx dfld'^q{n/qf is then the po- 
tential energy of a single particle- hole pair excited above 
the ground state. In this language, an increase of 
(5X2 {T, H) with both H and T can be understood as the 
consequence of the fact that the magnetic field gaps out 
soft particle-hole pairs, suppressing their contribution to 
the thermodynamic potential. 



FIG. 3: Skeleton backscattering diagram. 




B. Beyond second order 

Higher-order diagrams for S(T, H) can be divided into 
two groups. The first group is formed by diagrams in 
which a nonanalyticity is produced in the same way as 
at second order: by extracting a product of only two dy- 
namic bubble from the whole diagram. The rest of the 
diagram goes into dressing up of the fermion propaga- 
tors and renormalization of the 2fci? interaction lines into 
full static vertices. In real-frequency language, these di- 
agrams describe higher-order corrections to the effective 
static interaction in a single-pair processJ^ The second 
group is formed by diagrams in which a nonanalyticity is 
produced by combining more than two dynamic bubbles. 

These two groups of diagrams describe two distinct 
physical processes. As we will show in this Section, 
the first group corresponds to scattering events in which 
fermions, moving in almost opposite directions before the 
collision, reverse their respective directions of motion. 
We dub this process as " backscattering" . The second 
group describes scattering events with no correlation be- 
tween initial directions of motion. 

Third-order diagrams 6, c, and / in Fig. [2] belong to the 
first group. Nonanalytic contributions to x from these di- 
agrams are obtained by selecting two dynamic up-down 
bubbles and setting q = ri„i = in the rest of the dia- 
gram. As an example, we consider diagram h. Fermions 
from any of the two bubbles with opposite spins can be 
re-grouped into two up-down bubbles in the same way 
as in the second-order diagram in Fig. [TJ Retaining only 
the dynamic part of these two bubbles, we obtain the 
same nonanalyticity as at second order. The remain- 
ing, third bubble can then be evaluated at zero external 
frequency, which means that it renormalizes the static 
vertex. Diagrams of the first type to all orders can be 
cast into a single skeleton diagram, shown in Fig. [3l The 
fermion Green's functions in this diagram are of a FL 
form 



where is the renormalized Fermi velocity, Z is the 
quasiparticle residue, and A = 2(1^11 with jig being 
the effective Bohr magneton, which we discuss below. 
A hatched block in Fig. [3] is the spin component of the 
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renormalized static vertex, rs(k, p; q), obtained from the 
dynamic one in the hmit of Qm/^pQ ~* 0. We will follow 
a standard procedure^ and absorb factors of Z into Tg. 
In the low-energy limit (T, A <C ep), the fermion mo- 
menta k = rifcfc and p = ripp are confined to the Fermi 
surface, so that Ts depends on the angle between and 
Up-, (as well as on q) : Tg = Tg (n^ ■ np;q). At first order 
in U, the vertex reduces to Tsik, p; q) = — ?7(|k— p|), and 
only U{2kp) contributes to the nonanalyticity. Beyond 
the lowest order, however, more complicated angular av- 
erages of the interaction occur, and it is not a'priori clear 
what the prefactor of the nonanalytic term is. We now 



show, using the same procedure as in Ref. [ ^T§\ , that this 
prefactor is precisely the square of the spin component of 
the backscattering amplitude: i^'^TKuk • % = — 1, g = 0). 



1. Contribution to the susceptibility from diagrams with 
two dynamic bubbles 



The nonanalytic contribution of the skeleton diagram 
in Fig. [3] is given byi^ 



^2s(T,i?) = -^T^ J duk J dnpX[iyTg (uk ■ np;q)f 'Pu{^m,q;nk)'Pu i^ni,q;np), 

I 



(15) 



where 



'Pii {^m,q; nfc) 



ifljn/vpq 



2tt 



iil,n + A ) /vpq - rifc • 



(16) 



is the propagator of a particle-hole pair moving with in 
the direction of Uk with small energy f2m and momentum 
q. [P-\ [ in Eq. lU is obtained from V^i by averaging over 
rife: P||(17™,g) = JdnfeT'tx (f7,„,g;nfe).] The vertex F^ 
can be expanded in angular harmonics as 

oo 

Tg (nfe • rip; (j) = ^ F^,; {q) cos (Wkp) , (17) 

where 9kp = cos~^ (rifeHp). Substituting this expansion 
into (dni), we obtain 

2 ^ \vpqj 
xJ2^sM^sj,{q)Aw, (18) 

1,1' 



where 
Aw 

and 



(27r) 



J j dnkdnpCos{Wkp)coa{l'ekp) 



— rife • Uq IZ — Up ■ Uq 



rirn. - iA 



vpq 



(19) 



(20) 



Using the identities 

{ia-by^ = -zsgn(Rea) / dAe-^"sn(Rca)(a-<b) 

Jo 

i^JiiX) = I —e'^'°'\osiW) (21) 
Jo 



and 



dXJi (A) e 



-aX 



1 



(valid for Rea > 0), we obtain for Aw 



(22) 



A,, 



(-) 



i+i' 



(x/i^-z)'"' +(Vi^-z) 



Z2 + 1 



The expression for (T, H) reduces then to 
S2.(r,ff) = -iT^i?(z) 



- ^A)2 + {ipqf 



(23) 



(24) 



where 



B{z) 



E 



y.{-)'+''v^Tgj{q)Tgw{q). 



(25) 



As it was the case for the second-order diagram, the 
nonanalyticity in {T, H) is associated with the log- 
arithmic divergence of the integral over q [see Eq. ([8])]. 
Because the logarithm comes from the "tails" of the in- 
tegrand, typical Vpq are much lager than both and 
A, i.e., typical z are small. Therefore, one can safely put 
z = in the factor in square brackets of Eq. , upon 
which it reduces to unity. On the other hand, since q is 
still smaller than the momentum cutoff of the interaction, 
one can set q = in the vertex. Next, we recall that the 
small- momentum limit of z/Fs(k, p;q) is the scattering 
amplitude fs{dkp) ( RefQ). Therefore, 



B{0) 



.1=0 



(26) 
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which is a square of the exact backscattering amphtude. 

The rest of the integral in Eq. (|25p is evaluated in the 
same way as it was done at second order and yields the 
same scaling form as in Eq. pT|) . with vp ^ vp and A 
A. Therefore, the contribution of the skeleton diagram 
in Fig. [3] to the spin susceptibility is given by 



SX2S {T,H)=[U7t)Y 




(27) 



where ep = vpkp/2 is the renormalized Fermi energy. 

The renormalized Fermi velocity vp and Bohr mag- 
neton fls can be expressed in terms of the Landau pa- 
rameters gd and gs^i (Ref. 1), where c and s stand 
for charge and spin. The Fermi velocity is given by 
Vp = vp/{l + gc,i), while renormalization of the Bohr 
magneton follows from the requirement that the Zee- 
man energy of a spin in the magnetic field is 2flBH = 
2fiBH/{l + gs,n)- Then, 



fJ-B 



fJ'B 



1 + 5. 



s,0 



(28) 



We also recall that harmonics of the Landau interac- 
tion function are related to harmonics of the scattering 
amplitude^. In 2D, this relation is given by 



9a,' 



1 + 



ga..i 



(29) 



2-Sn.O 



where a = c, s. To first order in the interaction. 



2[/(0)<5„,o 
d0 



d6 

— cos nOU {2kp sin6'/2) 



Expressing via x(0, 0) and substituting the result 
back into Eq. ((27|) . we obtain 



(32) 



Before concluding this Section, we note that the ex- 
act backscattering amplitude /s(7r) depends logarithmi- 
cally on T and H due to singular renormalizations in the 
Cooper channel To see this, one needs to re- 
call that fs (tt) is equal (up to a prefactor) to the Cooper 
vertex for scattering from the states with momenta k and 
— k into the states with momenta — k and k, respectively. 
We will discuss this special feature of the backscattering 
amplitude in Sees. Ill C II and III C 21 but for a moment 
continue with the consideration of higher-order contribu- 
tions to s(r, H). 



Contributions to susceptibility from more than three 
dynamic bubbles 



cos nOU {2k p sin 9/ 2) . 



There are other diagrams at third and higher orders, 
which do not belong to the skeleton diagram in Fig. [H In 
zero magnetic field, these additional diagrams__;^ield only 
analytic contributions to S(T,iJ = 0) (Refs. S^^Ts'.'sS). 
This is not so in the presence of the magnetic field, as we 
are now going to demonstrate. 

At third order, there is only one diagram which cannot 
be fully absorbed into Fig. ([31)-diagram e in Fig. For 
a local interaction {U{q) = const), this diagram contains 
(30) a cube of the up-down bubble 



In general, g 
Z^{vplvp)vTli, 



where T" and r« are renormalized 
vertices at small momentum and frequency transfers, in 
the limits vpg/ftm and flm/vpq 0, respectively. 

It is convenient to re-express 6xiT, H) in terms of the 
actual (renormalized) spin susceptibility at T = _ff = 0, 
rather than of the susceptibility of a Fermi gas Xo^ ■ The 
renormalized spin susceptibility at = T = is given 



X(0,0) = 



m Vp 



Ms 



Vp 



^ -, , = — (1 - fs,o)Xo 

TT vp 1+ gsfi vp 



2D 



(31) 



S3e(r,i7) = --T^(i-Pu) 



(33) 



As one can readily verify, the first two terms do not give 
rise to non-analyticities, while the P-^^ term has already 
been accounted for in the skeleton diagram of Fig. 
The new contribution comes from the Pj^ term. Keeping 
only this term and integrating over q, we obtain 



E3e{T,H)^—TY, dqq- 



{n^-iAy + vW 



3/2 siTvy ^^^ni + A^' 



(34) 



n,„>o 



Subtracting off the ultraviolet contribution and summing 



over rJm, we find 



S3e(T,i?) = 



FT A h^coth — . 

Unvy ' 2T 



(35) 
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Differentiating twice with respect to the field, we obtain 
the new contribution to the susceptibihty 

5X3e(T,i/) = -u3^i?(^^)xr, (36) 

where 

2 

R{x) = coth X : — 2 \ r— 3— • (37) 

sinh X 3 sinh x 

In the two hmits, R{x — > 00) = 1 and R{x — > 0) = 1/Sx. 
We see that 6x3e has the same nonanalytic dependence 
on T and H as the second-order diagram: it scales lin- 
early with the largest of the two energy scales 

3 max{|A|,2r/3} 

0X3e(T,H) = -u . (38) 

eF 

There is one essential difference between the second- 
and third-order contributions: the nonanalyticity in 
6x3e (T, H) does not arise from a logarithmically diver- 
gent integral over q. Indeed, the momentum integral 



in Eq. ((34)) is convergent and comes from the region 
q ~ \^m\/vF ~ |A| /vf- This means that Eq. ([36]) can- 
not be obtained by replacing the dynamic part of the 
bubble by its asymptotic form at large vpq/flm, which 
was the case for the backscattering contribution. 

Notice that the sign of the third-order non- 
backscattering contribution is opposite to the second or- 
der result. This opens a possibility of inverting the sign 
of 6x in the non-perturbative regime (see Sec. Ill Cl for a 
more detailed discussion). 

To go beyond the perturbation theory for this new 
type of processes, we apply the same procedure as for 
backscattering. Namely, we combine all diagrams with 
three dynamic bubbles into a "third-order" skeleton dia- 
gram by replacing the bare interactions in Fig. [Sje by the 
renormalized vertices evaluated in the limit of flm/vpQ ^ 
0. This limit ensures that we obtain contributions with 
no more than three dynamic bubbles. The renormalized 
vertices are then again the spin components of the scat- 
tering amplitude fs , and the third-order skeleton diagram 
reduces to 



^3s,3 



-^X! y"^"fe J dup J dnsfs{nk-ni)fs{ni-np)fs{ni-np)Vu{^m,q;nk)Vui^ni,q;ni)Vni^ni,q]np){39) 



As it was done for the second-order skeleton diagram, 
we replace the bare vp and A by their renormalized val- 
ues and absorb the quasiparticle residue Z into fs- We 
now show that the three vertices in Eq. (|39p do not form 
the cube of the backscattering amplitude, i.e., that the 
nonanalyticity in Sxss comes from scattering of fermions 
with uncorrelated directions of the initial momenta. To 
demonstrate this, we adopt a simplified model, in which 
the angular dependence of fs (n; - rip) is approximated by 
the first two harmonics: 

fs{nk ■ lip) = fsfi + Uk ■ Upfs^i. (40) 



In this model, the backscattering amplitude is equal to 



/s(^)=/s.O-/s,l=/s,o(l-^)- (41) 



Substituting Eq. into Eq. performing straight- 
forward angular integrations, and differentiating twice 
with respect to the magnetic field, we obtain for the 
three-bubble contribution to the spin susceptibility: 



\Js,o/ \fJ-B/ ep 



5x3s{H = 0,T) = {UofF,{^ 

Js,0 



fJ-B J £F 
2 



AO 



Ms 



2D 



(42a) 
(42b) 



where Obviously, the product (/s,o)'^^3 (/s,i//s,o) in the pref- 

^ , , , „ r, -,\ „ 2 actors of Eqs. (I42al42bl) does not reduce to the cube of 

Fs{x) = l-3(21n2-l)a. + 3(31n2-2)x /,(^) from Eq. gj). 

+ (5/2 -3 In 2) (43) 
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A similar consideration can be extended to higher or- 
ders. At fourth order, we get an additional nonana- 
lytic contribution to x from processes with four dynamic 
particle- hole bubbles, at fifth order-from five dynamic 
bubbles, and so on. Each of these contributions can 
be converted into a skeleton diagram by dressing up 



fermion Green's functions and interaction lines, and nei- 
ther of them is expressed solely via the backscattering 
amplitude. For example, approximating /s(nfe • rip) as in 
Eq. ([301), we obtain for the fourth-order skeleton contri- 
bution 



\js,0 



<5X4.(i/ = 0,r)-(/..o)'F4(fi 

/s,0 



(44a) 
(44b) 



where 



1 - 4(3 - 4 In 2)x - (50 - 6 In 2)x^ 

41 



'(50 - 721n2).x^ 



20 In 2 



x4.(45) 



3. Isotropic scattering 

If one further neglects fs.i compared to fsfi, i.e., ap- 
proximates the scattering amplitude by a constant, the 
contributions to the thermodynamic potential from skele- 
ton diagrams from all orders form geometric series and 
can be summed up. Doing so, we obtain 



E{T,H 



lAI 



24:77 {vf)'^ 



2471(5^^)2 Vl-/s,o 
for |A| > T, and 

TA^ /I 



24:Tr{vF) 



••) 

(46) 



TA^ 



47r(wF)^ 

rA2 



ln(l-/,,o) + 



/s,0 



47r(wF)^ 
for lAI < T. 



ln(l + .g5,o) ^ + 9sfi 



(47) 



Differentiating ^(T, H) with respect to the field, we 
obtain 



6x{T = 0,H) 



and 



2 (t^\\M 2D 



5x{T,H^Q) 



ln(l-/,,o) + 



fsfi 



s,0 



(48) 



As \ ^ , 2D 



ln(l + 5s,o) ^+5s,o 



As \ T 2D 



(49) 



Equation (|49l) , without FL renormalization of the Fermi 
energy, was derived earlier in Refs. [TgHTtIITsI . 

In what follows, we will also need a full expression for 
the spin susceptibility for the case in which the angular 
dependence of the scattering amplitude is approximated 
by first two harmonics, as in Eq. (j40p . Such an expression 
can be obtained for the case when f^ i <^ f^^o while fs^o 
is arbitrary. The calculation of Sx is tedious but straight- 
forward. We present the result only for 6x{T, H = 0) 



SxiT.H^Q) 
where 



2^i^i(/s.o) 

Js,0 



Ms 



(50) 



Fo{x) 
Fi{x) 



In (1 - a;) + 



1 - X 



41n2 41n(l-x) 2 2{x'^ + 1)' 

(a; -1-1)3 ~ {x + 1)^ ~ (a; + 1)2(1 -x) " {x^ - 1)2 _ 



(51) 
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In the two limits, Fi{x < 1) « and Fi(x > 1) « 2x. 



C. The sign of the temperature and magnetic- field 
dependences of the spin susceptibility 

In the previous Section, we calculated the third and 
fourth order skeleton diagrams for a model form of fs {9) 
given by Eq. (|40)) . Beyond weak coupling, the expansion 
in skeleton diagrams does not have a natural small pa- 

I 



with functions F3{x) and ^4(0;) defined in Eqs. (|43l) and 
(US)) , respectively. The first term in Eq. ([55)1 is the square 
of the backscattering amplitude in the two-harmonic ap- 
proximation [cf. Eq. dH])]. 

Likewise, the T dependence of Sx{T, H ~ 0) is given 
by the sum of Eqs. (P7|) [taken in the limit oi H ^ 0], 

I 



In Figs. [3] and O we plot S'//(/s,o, /s,i) and 
Srifsfl, fs.i), correspondingly, as functions of fs,i/fs,a 
for a range of /s^o- We see that if both |/s.o| and |/s.i//s,ol 
are sufficiently large (but still less than one), the signs 
of the slopes are opposite to those of the backscattering 
contribution, i.e., x decreases with T and \H\. 

To get an idea about the numerical values of fsfi 
and fs^i, we use available data for Landau parame- 
ters. A system of fermions with repulsive interaction 
is expected to exhibit enhanced ferromagnetic fluctua- 
tions, which corresponds to a negative value of g^.o- In- 
deed, the Landau parameter g^^Q is negative in He'^ (in 
both bulk^ and film^^ forms), 2D gases in semiconduc- 
tor heterostructureS j^^'^^i^^ and many other fermion sys- 
tems. In bulk HeS, gs,o = —0.70 and gss — —0.55 at 
ambient pressure. '^^ In Si MOSFETs, 53,0 is also close to 
—0.7 in a wide interval of densities 1^ The first harmonic 
of the spin Landau function, gs.i has not been measured 
in 2D gases. Taking the bulk He'^ values as rough esti- 
mates for the 2D case as well, we obtain with the help of 
Eq. m- fs,o = -2.3, = -0.76, and fs,i/fs,o = 0.33. 
Although the magnitude of fs^ is probably too large for 
our truncated perturbation theory to be accurate. Figs. [3] 
and [S] indicate that both Sx{T,0) and Sx{0,H) are al- 



rameter. Still, it is worthwhile to analyze the result for 
not too strong interaction. To fourth order in the scat- 
tering amplitude, the field dependence of Sx{T = 0, H) 
is given by the sum of Eqs. ([77|) [taken in the limit of 
T ^ 0], (|j2a| . and (pa|) . Explicitly, 

xiT = 0, H) = ^ (^) ' x^Sh if so, fs,i) , (52) 
where 



(53) 

I 

(|42bl) . and (|44a|) : 

X(T, H = 0) = ^ (^] ' x^St if so, fs,i) , (54) 
where 



(55) 

I 

ready negative for this value of fs,i/ fs,o- 

We thus see that the field and temperature depen- 
dences of X are non-universal: while the slopes are posi- 
tive at weak coupling, they well may become negative at 
sufficiently strong coupling. [Later on, however, we will 
show that in the vicinity of a ferromagnetic QCP the sign 
of SxiT, H) is definitely positive.] 



1. Cooper renormalization 

References [TgHTsI considered a more subtle mechanism 
for changing the sign of 5x as compared to the second- 
order result, namely, renormalization of the backscatter- 
ing amplitude in the Cooper channel. As we have already 
said, the backscattering amplitude is special in that it 
is given by a fully renormalized vertex with zero total 
incoming momentum and momentum transfer of 2kF- 
Therefore, it can be expressed via angular harmonics of 
the irreducible Cooper amplitude, 7c, as 

^ l + "fc,i^T^iW/E) 



Sh (/sO, /s,l) — 77 ifs,0 - fs,iy + f 3,0^3 ( 

Js,0 



f4 TTi / fs,l 
nfsfi^i 1 

/s,0 



St (/sO, /s,l) — 77 (/s,0 - fs,lf + 



fsfi 



-JloF. ( 

Js,0 
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where E = max{T, A} is an appropriate energy scale. 
This gives rise to two effects. First, if at least one 
of 7c,i is negative, i.e., jc.io < 0? the system under- 
goes a superconducting transition of the Kohn-Luttinger 
type'^'^ into a state with orbital momentum /q at E'kl — 
M^exp (— |7(7,io|). The backscattering contribution to 
the spin susceptibihty 5x2s oc E/ {1 + '-fcdo^n{W/E)) 
diverges at Ekl as well. Above Ekl, Sx2s is non- 
monotonic: it decreases with E for Ekl < E < e^E Kj. y 
7.39Ekl and increases with E ioi E > g^Ekl (Ref.llGl). 
For E » e^ExL, Cooper renormahzation is weak and 



i5x2s becomes linear in E. Second, if all jc.i are positive, 
/s(7r) scales down to zero as 1/ ln{W/E) for — > 0. Con- 
sequently, the backscattering contribution to x is reduced 
by a factor of l/\n^{W/E). In this situation, the domi- 
nant contribution to Sx comes from non-backscattering 
terms, which do not contain singular Cooper renormal- 
izations. In Ref. [l^ this effect was accounted for in a 
model of isotropic scattering amplitude, fsid) — fs.o, 
by subtracting off the backscattering contribution from 
Eqs. dm) and (gH]). This gives 



SxiT = 0,H^O) 



fs,0 



fs.Q 



|A| 2D 



f^B \ 1^1 2D 



i^xr (57) 



and 



Sx{T^O,H = 0) 



ln(l-/. 



fs,a _ 1 f2 
l-/,s,o 2^^^o 



/s,0 



Mb\ T 2D 

— — Xo • 



(58) 



The signs of the H and T dependences in Eqs. ([57]) and 
(|58p now coincide with the sign of /s.o! for negative fs^, 
expected for a repulsive interaction, they are opposite to 
the second-order result. This mechanism was proposed 
in Ref. [H as an explanation of the negative sign of the 
slope of Sx{T, H = 0)(xT observed in Ref. HI 

A simple way to estimate the validity of the approxi- 
mation used in Eqs. ([57| and ([58|) is to consider the scat- 
tering amplitude with two rather than one components: 



fs{(f>) ^ fs.O + fs.lCOS(f>. 



(59) 



The vanishing of /^(Tr) at T — > implies that f^fl = fs,i 
at T = 0. In the approximation used to derive (jSSp . 
this relation was accounted for in the quadratic but not 
in higher-order terms. Substituting /s.o = /s.i into the 
third and fourth order terms in /g, we obtain, instead of 
Eq. dill). 



(0.02) X |/3o + (-45.65) x -J^^ + 



5xiT,H = Q) = 



in ([M)) are negative for fs^ < 0; hence, to this order, 
Sx{T, H = 0) < 0, which is consistent with Ref. 18. 



2. Coulomb interaction in the large N limit 



(60) 



Comparing ((60|) to ([58]) . we see that the prefactors dif- 
fer substantially, making it difficult to draw a general 
conclusion. At the same time, the signs of both terms 



Another issue is that the r esults of Refs.[il[il, as well 
as Eq. ([SO)) , are valid only below a characteristic en- 
ergy scale at which Cooper renormalizations of fsfi be- 
come significant. For a weak interaction, this scale is 
exponentially small. To estimate this scale beyond the 
weak-coupling regime, we consider the effect of Cooper 
renormahzation of the backscattering amplitude on the 
spin susceptibility in a large N model for the Coulomb 
interaction, developed earlier in Refs. l40H4lll4^ . To be 
specific, we assume that there are iV„ degenerate elec- 
tron valleys, so that the total (spinxvalley) degeneracy 
is A'' = 2A"„. While this model is especially relevant to 
Si- and AlAs-based heterostructures, which have at least 
two valleys {N = 4), it can also provide a useful insight 
even for single- valley system {N = 2). 

In an A^-fold degenerate 2D Fermi gas, the Fermi 
momentum is scaled down by a factor of ^/N : kp = 
ji/^/IttN, where n is the number density of electrons. 
On the other hand, the inverse screening radius k, which 
is proportional to the density of states, is scaled up by a 
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FIG. 4: (color online) Function Sh ( fs,o, .fs,i), which deter- 
mines the sign of the magnetic-field dependence of the spin 
susceptibility calculated to fourth order in the skeleton in- 
teraction [cf. Eq. (|53|) ]. plotted as a function of /s,i//s,o 
for fs,o = -0.3 (dashed), -0.5 (dotted), -0.7 (dashed)'. The 
curve for fs,o — —0.3 was multiplied by a factor of 10 for 
clarity. 



is only a weak-coupling regime. If A'' ^ 1, there are two 
regimes: weak coupling (g^v *C 1) and strong coupling 
{gN 3> 1); the latter is of the most interest for us. 

The details of the calculation are given in Appendix [Xl 
Here we present only the result for the field-dependent 
part of the spin susceptibility 



o 



X -3 



CO 




6x 



Sep \L, 



N 



AO I 



(62) 



where Lc = ln(ei?/|A|) is the Cooper logarithm and, as 
before, A = 2^bH. Equation ([^ is valid for Lc ^ 
1, i.e., for A <C Ec = e_Fexp(— A^). The first term in 
Eq. ([S^ is the backscattering contribution, which is the 
leading term in the \/N expansion. The second term 
is the contribution from other processes, which is the 
next-to-leading term in this expansion. As expected, the 
backscattering contribution is scaled down by a factor of 
1/L^. The change of sign occurs at Lc = N^^"^ /\/2 or 



A = £;* = cxp (-TV^/V^) < Ec- (63) 



FIG. 5: (color online) Function St {fs,o, fs,\), which deter- 
mines the sign of the magnetic-field dependence of the spin 
susceptibility calculated to fourth order in the skeleton in- 
teraction [cf. Eq. (|55p ]. plotted as a function of /s,i//s,o for 
fs,o = -0.3 (dashed), -0.5 (dotted), -0.7 (solid). The curve 
for fsfi — —0.3 was multiplied by a factor of 10 for clarity. 



0.5 




^ -0.2 -0.1 0.1 0.2 0.3 



f 7f 

S,1 S,i 



factor of N. Their ratio, 

gN = n/kp = r,iV3/2/2 (gi) 

with rg — me'^ / y/iTn, defines an effective coupling con- 
stant. We still need to assume that ^ 1; only then the 
mean-field, random-phase approximation (RPA) is valid. 
If iV ~ 1, then gN ^ Ts 1, which implies that there 



For an estimate, we consider a 2D electron gas in the 
(001) plane of a Si MOSFET, where = 4; correspond- 
ingly, Ec = OMSep and E* = 0.0035eF. In the experi- 
ment of Ref-IU, the highest Fermi energy in this measure- 
ment is about 40K. Then, Ec = 0.70 K and E* = 0.14 
K. Both energies are smaller than the disorder broaden- 
ing in these samples. This shows that the mechanism 
of the sign reversal of 6x due to Cooper renormalization 
does not have room to develop until disorder becomes 
important. Sx in Ref. [l^. 

Notice also that Ec is still larger than the energy 
scale of the Kohn-Luttinger superconducting instabil- 
ity. Indeed, in 2D the Kohn-Luttinger effect starts only 
at third order in the interaction,"^^ which implies that 
Ekl exp(-iV3)eF < Ec- 

The above estimates are based on additional assump- 
tions, such as ^ 1, and therefore cannot give rigorous 
results regarding the T and H dependences measured in 
2D heterostructures. The crucial experimental check for 
the many-body nature of these effects is the T/H scaling 
of the susceptibility which, to the best of our knowledge, 
has not been performed in detail. 

We should also point out that while the T-dependence 
of X in Si MOSFET was obtained in a thermodynamic 
measurement (via the magnetocapacitance) the H- 
dependence of x is all heterostructure a^^'^^'^^ was ex- 
tracted from Shubnikov-de Haas oscillations. While 
it is known that Shubnikov-de Haas oscillations con- 
tain a renormalized spin susceptibility of a FL, x(0, 0), 
it remains to be verified that the field-dependence of 
the susceptibility can also be extracted from such a 
measurement 4^ 
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D. Thermodynamic potential of a 2D Fermi liquid 
as a function of magnetization 

Preparing the ground for the analysis of a ferromag- 
netic QCP in Sec, mil it is convenient to obtain the ther- 
modynamic potential in terms of magnetization rather 
than of magnetic field. In this formulation, the suscepti- 
bility is defined as 



X 



1 a^s 



(64) 



is the number density of 



where M — n-^ — and n-^ 
spin-up/down fcrmions. 

In the RPA, which neglects FL renormalizations {Z — 
vp/vp — 1), the recipe for finding the free energy was 
given in Ref. [s^. For a local interaction U = u/v, 



E{A,M,T) = - 



vA^ MA uAP 



A 2 Av 

where (5S(r, A) is the sum of the RPA (ladder) series 



(65) 



<5s = T V (In [1 + (V^)nn] - («/^)nTx) , 
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(66) 



(The second term in Eq. (|66p compensates for the first- 
order contribution not present in the ladder series.) The 
relation between M and A is found from the condition 
(dE/dA) 1^^^ = 0, which gives M — vA. Neglecting the 
RPA term, one obtains the Stoner-like spin susceptibil- 
ity x(0, 0) — 2^\v/{l — u), which is consistent with |3T|) 
for Z = vp/vp = 1 and /o,^ = —u/{l — u) (or, equiva- 
lently, — ~u). Evaluating further x(r. A) with the 
RPA term included, and using the relation between A 
and M, one reproduces Eqs. pS)) and (^5]) without FL 
renormalizations. 



Fermi-liquid renormalizations 



Since our primary interest is the spin susceptibility at 
small momenta, we focus on the Pomeranchuk instabil- 
ity towards a ferromagnetic state. In the FL theory, this 
instability occurs when gsfl approaches —1. All other 
partial components of the spin and charge scattering am- 
plitudes are assumed to remain finite at criticality and, 
without loss of generality, can be taken to be small. 

The primary goal of the present Section is to demon- 
strate that some terms in the thermodynamic poten- 
tial of Eq. (|65p are renormalized at energies compa- 
rable to the bandwidth, where S is static, while oth- 
ers are renormalized at much smaller energies, of order 
Ta.SiK{T ^ M / [v pv / V p)} , where is dynamic. 

Consider first the A^ term in Eq. ([65]) . which is the 
thermodynamic potential of free fermions. For a FL, this 
term can be calculated using the renormalized Green's 
function Applying the Luttinger-Ward formula^ 

for the thermodynamic potential of free fermions and ex- 
panding S to order A^, we obtain at T = 



:(0) _ _ 



E 

vA^ 



doJrn 



w 

-w 
w 



vpd{k~ kp)\n{-G-^ 
vpd{k — kp) 



-w {vpi^k — kp) — icOm) 



(69) 



The integrals in Eq. (I69p are controlled by large momenta 
and energies, of order uj„i ~ vp\k — kp\ ~ W. Integrating 



first over k and then over cjm, we find that 



Vp vi^ 

vp 4 



(70) 



To evaluate the second term in ((65|) . we need the re- 
lation between the magnetization and Zeeman energy, 
which can be found by expressing the number densities 
of spin-up and spin-down fcrmions in terms of the Green's 
functions 



Equation can be generalized to a FL. In this Sec- 
tion, we assume that the system is away from the im- 
mediate vicinity of a QCP, and the effective interaction 
can be considered as static. We discuss specific condi- 
tions below. For a static interaction the relevant fermion 
self-energy depends on k but not on ujm'- 



Yf{uJm, fc) « - {vp - Vp) (fc - kp). 



(67) 



In this case, the quasiparticle residue Z = 
{1 — idY, / duj„i) ^ is equal to unity and the fermion 
Green's function is given by 



GT,i{ujm,k) 



1 



vp{k - kp) ± A/2 



(68) 



where, as before, A = 2flBH — 2fiBH/{l + gs,o)- We 
will see later in Sec. lIIII that the self-energy becomes pre- 
dominantly Lo^- but not k- dependent in the immediate 
vicinity of a QCP. 



M = nt 



~2^ 



vpd{k~kp){G^~Gx\. (71) 



Using the Green's functions from Eq. , we again find 
that the integrals in Eq. (|7ip are controlled by energies 
of order W . Performing the integrations, we obtain 



M = .^A. 

Vp 



(72) 



Using ([70]) and recalling that the relation (17^ must follow 
from the condition 9S/9A = 0, we find that the second 
term in (j65p retains its form, but A changes to A. 

In a similar way, the Hubbard U in the Hartree term 
in Eq. (j65p is replaced by the Landau parameter 



U ■ 



u ^ _vf_ gs,o 

V Vp v ' 



(73) 
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so that the Hartree term becomes 



parameters gsfl and vp/vp- At T = 0, 



As a result, we have 



B(T,M,A) = 



VF 



2 AvFi- 



'V 



M'^. (74) 

l-<5S(T,A). 
(75) 

So far, aU renormahzations in Eq. (|75p are from energies 
of order W . 

Consider next the RPA term, (5S(T, A). Re-calculating 
it with G{k, ojm) from Eq. we find that it retains the 
same form as in Eq. ((66l) , except that U is replaced again 
by the Landau parameter and the polarization bubble 
lit I contains the renormalized Fermi velocity: 



SE{T, A) = T 



9 



In 



9s,qVf 
vfv 



n 



Ti 



9sflVF 

vfv 



n 



Ti 



There are both analytic and nonanalytic terms in 5^. 
The analytic, A^ contribution comes from energies 0{W) 
and can be absorbed into the NP term in Eq. (|75p . once 
the prefactor is expressed in terms of the Landau param- 
eter Qsfi rather than of the bare interaction. The leading 
nonanalytic term has the same |Ap form, as in the per- 
turbation theory, but the prefactor now contains the FL 



5S(0,A) 



s,0 



(77) 



The key point is that the A^ term in Eq. (|77p comes from 
small energies: uj^ ^ VF{k — kp) ^ A <^ W . Therefore, 
vp in Eq. (j77p is the Fermi velocity on the small-energy 
scale. 

The equilibrium condition {d'E,/dA)\^ = now 

gives M = Aiy{vp/vp) + 0{A'^), consistent with Eq. ([7^ . 
and the spin susceptibility at T = Af = is obtained 
from Eq. 



n\ VP 1 2 vp 

X(0,0) = — — = fig— — (1- /s,o) • 

TT Vp I + gs TT Vp 

(78) 

As expected, this result coincides with the general ex- 
pression for the renormalized spin susceptibility in a FL 
(Ref.lH). 

Using Eqs. ((75|) and ((7^ . we can now construct an 
expansion of the thermodynamic potential in powers of 
magnetization. To order |Mp and at T = 0, we obtain 



S(r = 0,M) = , , ^ M^-—^^(-^-y\M\^ + bM^ + ... 

= ^^(1 + 5.o)M2 _ ai^v^ , ^ ^ 

WVp 2i'KV'^VpVp 

I 



For completeness, we added a regular Af^ term to :i. 

Evaluating x{T = 0, M) with the help of Eq. ([Ml) and 
using the relation between M and equilibrium A, we re- 
produce the linear in H term in the spin susceptibility, 
Eq. dm). 

We emphasize again that the Af^ and \M\^ terms in 
Eq. ([7^ come from different energy scales. The Af^ term 
comes from high energies of order W , and vp in this term 
is the renormalized Fermi velocity at energies of order 



W . The \M\^ term comes from fermions with energies of 
order M / {v pv / v p) = A, and vp in the \M\^ term is the 
Fermi velocity on that scale. 

At finite temperature, the expression for the thermo- 
dynamic potential becomes more involved. We present 
only the result and show the details of the derivation 
later, in Sec. pil D ip . where we compute S(r, Af ) in the 
spin-fermion model. To logarithmic accuracy, we obtain 



S(T,A//) 



ip 



vp 

_ A/2 Vp 

Av Vp 



fsfi 



1 - /.,0 



+ ln(l - fs,>) 



Vp T 
vp ep 



fsfi \ fip 
l-fsfi) hl&^Tu^vl \vp 



1 + 5s,o - (5s,o + ln(H- g^^o) ^) — — 

' Vp ep 



■ bAr + ... 

(80) 



r 



The key result here is that for T ^ M/ [vpv/vp) the ex- pansion of ^ in powers of M becomes analytic: the | Af|" 
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term is replaced by an analytic M'^ /T term. Simulta- 
neously, the prefactor for the quadratic term acquires a 
linear in T correction, which is just a nonanalytic tem- 
perature dependence of the spin susceptibility. Indeed, 
evaluating x{T^^), from ((80|) . we reproduce Eq. (|49)) . 



III. MAGNETIC RESPONSE NEAR A 2D 
FERROMAGNETIC QUANTUM CRITICAL 
POINT 

A. General considerations 

We now consider the immediate vicinity of a ferromag- 
netic QCP. In what follows, we first assume that a con- 
tinuous ferromagnetic transition does exists and obtain 
the thermodynamic potential S(M, T) along a continu- 
ous second-order transition line by extending Eqs. (|79|) 
and (|80p to energies below the scale where the self-energy 
crosses from static to dynamic forms. ( 46lf47[ ). Next, we 
show that this line becomes unstable at low enough tem- 
peratures because of nonanalyticities which survive even 
in the vicinity of the QCP. We will argue that the insta- 
bility may occur in two ways: i) the second-order phase 
transition into a uniform ferromagnetic phase becomes 
first order or ii) the transition occurs via an intermediate 
magnetic phase with a spiral magnetic order. More spe- 
cific predictions are possible within more specific models. 
One of such models is a model with a large radius of the 
interaction in the spin channel. ^■^ We will show that in 
this model the first-order instability occurs before the 
spiral one. 

A tendency towards the first-order transition can be 
seen already from Eq. ([79|) . Indeed, the cubic term in 
M in S(M, T = 0) is negative, which implies that a state 
with finite magnetization is energetically favorable. Close 
to the critical point, the term is small and the ther- 
modynamic potential (|79|) is negative over some range 
of M . This means that the first-order phase transition 
preempts the second-order one at T = 0. Indeed, for the 
thermodynamic potential of the form 



S(M) = aM^ - c\M\ 



(81) 



with c > 0, the magnetization jumps to finite value of 
Mq = 2a/c already for a — c^/Ab > 0, i.e., before the 
second-order transition takes place. 

At finite T, the |M|^ term in the thermodynamic po- 
tential crosses over into a —M'^/T one, which is still nega- 
tive. If T is low enough, this negative term is larger than 
the regular, term, and the transition remains first 
order until the —M'^/T term becomes smaller than the 
regular 6Af^ term. At higher T, the transition becomes 
second-order. 

This analysis is, however, incomplete because it is 
based on the result for S(T.M) derived under assump- 
tions that the quasiparticle residue Z = 1 and the effec- 
tive Fermi velocity vp — kp/m* is finite. As we have 
already mentioned, this is true only if the self-energy is 



static. Since the first-order jump in magnetization, Mq, 
is proportional to the critical parameter 1 -\- gsfi, the cor- 
responding energy scale Mq /v is also small and falls into 
the regime where the self-energy is dynamic and Eq. (j79p 
is no longer valid. 

If the self-energy is dynamic, the Z factor and vp/vp 
are both given by {1 — idT^{oJm) / duJm) (so that the 
product Zvp/vp remains intact.) This would not lead to 
substantial changes if Z remained finite at a QCP. How- 
ever, it is well established by now that d'S{uJm) / dcum di- 
verges at a ferromagnetic QCP in 2D; hence, both Z and 
vp/vp vanish.— i^i^i^iSSii^ One then might be tempted 
to conclude that the nonanalytic term in Eq. ([79|) van- 
ishes, as it is proportional as an overall factor to the 
renormalized velocity vp evaluated at low energies. We 
will show, however, that the nonanalytic term in the free 
energy survives even at the QCP, albeit in a weaker form 
(|M|3 is replaced by \M\'^/^). 

Before proceeding further, we mention two paradoxes 
with the vanishing of Z and vp/vp at a ferromagnetic 
QCP. First, there seems to be a contradiction with the 
Stoner criterion which says that a ferromagnetic transi- 
tion occurs at some critical, finite interaction strength. 
If we formally use the FL relation = Z{Zvp/vp)T^ 
(Ref. Q) with Z ^ Q but Zvpjvp = const, we find that 
the condition ggfi ~ —1 can be satisfied only if P^ oo. 
Second, in the FL theory, the velocity renormalization is 
determined by the I = 1 harmonic of the Landau function 
in the charge sector: v/vp = ifi/m — l-|-gc,i- Hence, the 
vanishing of vp implies that g^^i = oo. Meanwhile, the 
very idea of a Pomeranchuk instability is that it occurs 
only in one particular channel, e.g., in the spin channel 
with the angular momentum I = for a ferromagnetic 
QCP. All other channels, including the charge channel 
with / = 1, remain uncritical, which seems to be incon- 
sistent with the condition gc^i = oo. 

We make a few general remarks about these two para- 
doxes first and show specific results later. 

• The assumption of the conventional FL theory 
about a single relevant / = spin channel near a 
ferromagnetic QCP is valid if there is a wide range 
of energies below the cutoff W, where the fermion 
self-energy is static. Within this range, Z = 1 and 
vp/vp differs from its bare value only because of a 
non-singular interaction in the / = 1 channel. Then 
gs.o = Z"^ {v F / V f)vT^ is of the same order as i/P^^ 
and a critical value of gs.o = — 1 is approached al- 
ready at finite interaction strength. As we have 
already demonstrated, the Stoner enhancement of 
the spin susceptibility comes from fermions with 
energies of order W , hence, at energies below W , 
the susceptibility is already enhanced by the Stoner 
factor (1 -I- gs,Q)~^ > 1- 

• At some energy scale, A < W, the self-energy un- 
dergoes a crossover between static and dynamic 
forms. Accordingly, Z and vp/vp begin to vary 
below A and eventually flow to zero at the QCP. In 
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this regime, the conventional FL theory based on 
the static approximation is no longer valid and has 
to be replaced by a "new" low-energy FL theory, 
in which the "bare" fermions are the ones on the 
scale of A, the "bare" interaction between fermions 
is in the spin channel, and the interaction poten- 
tial U is replaced by the effective interaction, which 
scales as 1/(1 -I- (7s, o)- This low-energy FL theory is 
a spin-fermion model. The Landau parameters for 
the low-energy FL differ from those of the conven- 
tional FL. In particular, all harmonics gc^i, includ- 
ing (7c, 1 — vp/vp — 1 diverge at a 2D QCPi^ 

• The spin-fermion model is valid only if the crossover 
between static and dynamic forms of the self-energy 
occurs on a scale much smaller than W. Otherwise, 
one cannot consider only the I = spin channel. As 
we will see, the condition A W can be satisfied 
if the interaction is sufficiently long-range or, else, 
if the model is extended to ^ 1 fermion flavors. 
We will assume below that at least one of these two 
conditions is satisfied. 

Because both Z and vp/vp on the scale of A are just 
constants, we will absorb Z into the effective spin inter- 
action, and measure the velocity renormalization below 
the cutoff with respect to its value at the cutoff. In other 
words, we assume that the bare fermion propagator is 
G~^{k, ujm) = ii^m — vp{k — kp), and use symbols Z and 
vp to describe renormalizations at energies smaller than 
A. 



B. Spin-fermion model in zero magnetic field 

We now consider in detail the low-energy effective the- 
ory near a ferromagnetic QCP: the spin-fermion model 
We first review briefly the properties of this model in 
zero magnetic field^^ and then show how the model is 
modified in the presence of the field. 

The spin-fermion model includes low-energy fermions 
with a bare propagator G{k,uJm), collective spin excita- 
tions Sq, whose bare propagator is the static spin suscep- 
tibility x('?)i ^-iid the spin-fermion interaction, described 



by the Hamiltonian 



Hint — -TT 



9_ 
M 



k,q,Q:,a' 



(82) 



where A/" is the number of lattice sites. The spin-fermion 
coupling g is related to the Landau parameter g^fl as 
g — {TT/m)y/—gsfi. Near the QCP, gs.o ~ —1 and g « 
Tr/m. The bare boson propagator xil) is proportional to 
1/(1 + gs.o) for q ^ 0. We assume that the q dependence 
of X at small but finite q is described by the standard 
Ornstein-Zcrnike formula 



Xiq) 



1 



1 



where 



TT 1 + 5.,.o + (ag)2 ™2 ^-^ + q^' 



e ' = (1 + (7s,o)a~ 



(83) 



(84) 



Similar to the 1 -I- g^^ term, the analytic q^ term in 
xiq) comes from fermions with energies comparable to 
W. This term can be obtained in the RPA scheme, but 
one has to assume either that the dispersion is different 
from a free-fermion one, i.e., from k^/2m, or that the 
exchange interaction is momentum-dependent; otherwise 
the particle-hole polarization bubble does not depend on 
q for q < 2kp in 2D. If the q^ term is comes from the 
momentum dependence of the interaction, the length a 
is the radius of the interaction. In this case, the RPA is 
justified for a sufficiently long-ranged interaction, i.e., for 
a/cF > 1 (Ref.lHi.) 

The spin-fermion interaction affects both fermion and 
boson propagators. Collective spin excitations acquire a 
self-energy S^((7, il), while fermions acquire a self-energy 
S^(fc,ti;„j) which gives rise to renormalizations of Z and 
of the Fermi velocity: 



G{k,Uyn) 



x{q,^m) = 



1 



vp{k — kp) {k,L 



TT 1 +5^,0 + (09)2 + SB (g,l]„) 



.(85) 



To one-loop order, the T ~Q self-energies behave as 



vpq 

„F/, , X ( iX uJm, for w,„<Wo/A^, 

L (/C = A;F,Wm) = <^ 1/3, 2/3 r ,,/i3^,, .,,,1/4 3/4 



ILUq' U!,n , for UJq/X^ ^LUm ^ L^n 

^0 \ ^ 1 X 



^0 

i:'{k,uj^ = 0)^vp{k-kp)[j^j ^'^[j] (^6) 



where 



(87) 
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and S{x < 1) - X, S{x > 1) = 0(1). To simplify the 
formulas, we assume that uj„i > 0. Eq. ([86)1 is valid for 
t^o/A^ < Wmax, i.e., for (5 < 1. 

In the RPA, a ferromagnetic transition occurs at 5 = 1, 
but the critical value of g may differ from one in a more 
general model. For akp 1, luq is parametrically small 
compared to e^, i.e., the A;— dependent part of the self- 
energy is always smaller than vp{k — kp). 

To compare the frequency- and momentum- 
dependences of the self-energy, we consider the 
Green's function near the renormalized mass shell: 
vp{k - kp) = iujm + T.F{k,uJm)- For < uq, 

iojm ^ and vp\k — kp\ ^ uP'/^Luy^ near the mass 
shell. Consequently, the fc-dependent part of the self- 
energy is smaller than the cjm-dependent part. For 

U>rn > Wo, S^(Wm) < iuJra and tkVp\k - kp\ ^ LOjn 

near the mass shell. Comparing S^(fc,Lj„j = 0) and 
[k = kp,uj„i), we find that the two become compa- 
rable at U!rn ~ Wmax = ^^0 i'^p/^of^'^- For < W^ax, 

Y,^{k,u}m) depends predominantly on ojm while for 
ujjn 3> ti^max it dcpcnds predominantly on k. Note that 

2/3 

^max is also the upper cutoff for LUm scaling of the 
fermion self-energy (see Eq. (|86p ]. At larger ujm, 
scales as Inw^. 

The scale Wmax is larger than but still paramet- 
rically smaller than ep [indeed, Wmax ~ ujQ{akp)^ 
ep/{akp)]. The upper limit for the low-energy theory. A, 
can then be set somewhere in between Wmax and W ^ ep] 
its precise location being irrelevant as long as Wmax <SC ef. 

In what follows, we will also need the fermion self- 
energy at finite temperatures. At finite T, the fermions 
interact both with classical (fim — 0) and quantum 
{flm 7^ 0) spin fluctuations. The quantum contribution 
to the self-energy is a scaling function of uJm/T: 



where the scaling function Q{x) is such that Q{x ^ 1) = 
1 and Q{x ^ 1) ~ l/x^^^. The classical contribution 
contains a static propagator that diverges at QCP as 
q^^. This divergence can be regularized in two ways: 
by accounting for a thermally generated mass of spin 
fluctuations due to mode-mode coupling (not present 
in the spin-fermion model) ^^'^^ or by resumming the 
self-consistent Born seriesi^ In the flrst approach, the 
zero-temperature correlation length in Eq. ([55)1 ^(0) = 
a/y/l + gs,o is replaced by ^(T) - a/ [(T/A) In (A/T)]^/^ 
(Ref. [2l[ ). The classical part of the self-energy = 



parable at a characteristic temperature 



S {ujm = ttT, T) then becomes 



kpa? 



cx I 



T 



ln(A/T) 



1/2 



(89) 



Tqc = 9 [akp 



A3 



(90) 



For T <C Tqc the classical contribution dominates over 
the quantum one, and vice versa. 



1. Eliashberg theory 

We now focus on the low-energy region oJm ^ Wmax, 
where E^(fc, w„i) is predominantly dynamic. Within this 
region, there exists another scale, ojq ^ bj-^aax/ {akpY ^ 
'^max, at which the fermion self-energy becomes compa- 
rable to uj^. Below wq, 5]^(cjm) > Lo-^m i.e., the system is 
in a strong-coupling regime. Close enough to the QCP, 
i.e., for A 3> 1, the strong-coupling regime, on its turn, is 
divided into two more subrcgimes: i) loq/X^ < uJm < i^o, 
where the self-energy has a non-FL, uJm^, form and ii) 
ujrn < wq/A^, where the FL behavior is restored, i.e., 
ReE^(a;) oc u and ImE^(a;) cx w^. 

Since E^(a;m) > oJm at a;,„ < ujq, the accuracy of the 
one-loop approximation for the self-energy becomes an 
issue. Previous work""*'-^ demonstrated that the self- 
consistent one-loop approximation (the Eliashberg the- 
ory) cannot be controlled just by a large value of the 
parameter akp, as higher-order diagrams in the strong- 
coupling regime are of the same order in 1/akp as the 
one-loop diagram. To put the theory firmly under con- 
trol, one needs to extend it formally to 3> 1 fermion 
flavors; then higher order terms in the self-energy terms 
are small in {In N)/N'^. In what follows, we neglect this 
subtlety and assume that the Eliashberg theory is valid. 

The frequency-dependent self-energy T,^{uJm) from 
Eq. ([55]) leads to the renormalization of the Z-factor 
(equal to the inverse velocity renormalization factor). 
Right at the QCP, both Z and vp/vp depend on ujm 
as 



Vp 



1 - i 



duJrn 



= 1 



UJQ 



1/3N 



(91) 

The boson self-energy is generated by inserting the dy- 
namic fermion bubbles, made out full propagators, into 
the bare spin-fermion interaction. Summing up the RPA 
series for the renormalized spin-fermion interaction, we 
obtain 



(92a) 



(r!„-icoEF(17™))V4g2 



T7i92b) 



In the second approach, one obtains a self-consistent 
equation for E which yields a similar T dependence of 

The quantum and classical contributions become com- 



where is a slowly varying function of ^m, which inter- 
polates between two limits: cq = 1 for flm «C uJnJX^ and 
CO w 1.2 for LJo/X^ < f7„ < (Refs. Hilllzil) . For 
free fermions, P is the same as P introduced in Eq. ([5]). 
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In the limit of small frequencies, P{q, fim) reduces to 
the Landau damping form O^/wp-g. The static boson 
self-energy is small in l/(akp) and non-singular, and we 
neglect it. 

Equation l|92bp has to be treated with caution because 
P is the dynamic part of the particle-hole bubble made 
of dressed fermions but without vertex corrections. The 
latter are irrelevant for vpCl ^ ^ but are important for 
vpq ^, as they are necessary for the Ward identities 
to be satisfied. This problem is generic to all models in 
which the effective interaction is peaked at zero momen- 
tum transferi^ Fortunately, this complication does not 
arise in the study of non-analyticities in the thermody- 
namic potential and in the spin susceptibility because, as 
it will be shown later, we will only need to know P{q, ilm) 
for vpq ^ ^^m,S(r2m). Therefore, we will be using Eq. 
(j92bp in what follows. 

The thermodynamic potential of the spin-fermion 
model in the Eliashberg approximation was obtained in 
Ref.i (see also Sec HlTU]) : 



S(r) = SF {T) + -TJ2 



4^ 



\n[x-Hq,^m)] , (93) 



where Sp (T) = -vTY.^^^ \uj,n\ = -n^vT^ /'i is the T- 
dependent part of the thermodynamic potential of a free 
Fermi gas and x(9, ^^m) is given by Eq. ([85]) . Differenti- 
ating S(T) with respect to temperature, one obtains the 
specific heat C{T), which behaves as (1 -I- A)r away from 
the QCP and as T^/^ at the QCP. 



C. Spin-fermion model in a magnetic field 

We now return to our main discussion and consider the 
spin-fermion model in the presence of a magnetic field. 
First, we derive a general expression for the thermody- 
namic potential in a magnetic field, and then analyze 
the structure of the nonanalytic terms in the vicinity of 
a ferromagnetic QCP. For reasons already explained in 
Sec nil A] we take the bare fermion propagator as 



1 



VF{k- kF)± A/2 



(94) 



where A = and (Ib Mb/(1 + gs,o)- 



1. The Luttinger-Ward functional in a magnetic field 

We first derive the thermodynamic potential for the 
spin-fermion model in finite magnetic field, starting from 



the Luttinger-Ward functional^ and making use of the 
Eliashberg approximation, which neglects vertex correc- 

tion: 9 -57.58.59.60 

The Luttinger-Ward functional'*^ contains four terms 

^ S(M, A, T) = Sf + Sb + Sfb + Sm, (95) 
where .np is the potential of fermions dressed by the 
interaction with bosons, Sb is the potential of bosons 
dressed by the interaction with fermions, SpB is the skele- 
ton part which describes explicitly the fermion-boson 
interaction at low-energies, and 2m — gs^M"^ /{^v) + 
(l/2)AfA is an extra M-dependent high-energy contri- 
bution, same as in ([75]) . (There is no double counting, 
as one can verify explicitly.) For H = 0, the Eliash- 
berg form of the Luttinger-Ward functional was derived 
in Refs. [qI IstIIssIIboI . Extending the derivation to the case 
of finite spin polarization, we obtain 



TY,[\n{-G-'{k)/W}~J:l (fc)G.(fc)] 
^ (1 + 6,^,,) [in (g)} - EL' (q) D..' (q) 



T 
1 

9^ rr2 



Sfb = {^ + Sa,a')GAk)D^^,{k-k')G^,{k{m 



kk' 



where g is the dimensionless coupling constant defined 
in Eq. §7^, k = {k,ujm),q = (q,f^m), cr,cr' =T,i, 
and as before, summation over k and q implies summa- 
tion over Matsubara frequencies and integration over mo- 
menta. The functions G|^j are the exact fermion Green's 
function 

Gta = + ^li M - ek ± A/2) (97) 

and Dcra' (<?) is the propagator of spin fluctuations 

D-^^, {q) = l + gs,o + {aqf + S^,, (<?) , (98) 

where and E^^, are the fermion and boson self- 
energies, correspondingly. Notice that the fermion self- 
energy does not contain a constant part evaluated at 
ujm = and k = kp- this part has been absorbed into 
the renormalized Zeeman energy, A. 

By construction, the Luttinger-Ward functional is sta- 
tionary with respect to variations of the fermion and bo- 
son self-energies. The stationarity conditions 



6E 







(99) 



yield 
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9" 



EB^,(g) = ^Lr^G,(fc + 5)G,,(fc), (100a) 

A; 

~2 ~2 

(fc) = -^Tj2Da,a (q) {k + q)-^TY, [Da,-Al) + D-aAl)] {k + q) . (100b) 



In the presence of a magnetic field, there are two differ- 
ent boson self-energies: = and = E^^. E^^ 
is composed of fermions of the same spin (E^^ — E^^). 
It depends on the magnetization only via a shift of the 
chemical potential. This is a regular, analytic depen- 
dence which does not lead to nonanalyticities in the ther- 
modynamic potential. Neglecting this dependence, we 
set S'j^i = E^i = E^, where E^ is the boson self-energy 
in zero field, given by Eq. (|92ap . On the other hand, 
E|?j^ is composed of fermions with opposite spins and 
depends strongly on the magnetization via the Zeeman 



term: E?, = Ef|. = g^P-\i, where 



y (f^m - icoEF(r!„) - ia) +v%q^ 

(101) 

This also imphes that D-^^-^{q) — I?|_|(g) = D{q) and 
Dhlil) ~ P^l.lil)- Expressions for Sb and Sfb can be 
then simplified to 



Sb = r5][ln{i?-i(9)}-EB(Q)i?(g)]+|5][ln{i?f/(g)}-E^^(g)i?u('?) 



^FB 



2v 



kk' 



J^Dik- k') G„ (fc) G„ (k') + {k - k') Gt (fc) {k') 

(7 

I 



while the fermion self-energy changes to 



(102) 
(103) 



E^ {k) = -^^TY,D [q) G. {k^q)~ l.TY,Dn{<l)G-a {k + q) 



I 

The field-dependent part of the thermodynamic poten- tial involves only Dn (q) 
I 



(104) 



-T^ [In {-G;1 (fc) /W} - E^ (fc) G. (fc)] + I 5] [ini?-; (g) - E^ (q) Dn (q) 

k,a q 



-2 j\,,f2 1 

J2 Dn {k - fc') Gt (fc) Gx (fc') + + ^MA. 

kk' 



9 rp2 

2v 



(105) 



r 



With the help of Eq. (fTUDa|l . we simplify Eq. pUS)) to 
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M2 1 

4l7 + 2^ 



2 = +5.,o— + T^MA - [In {-G;^ (fc) /W^} - (fc) (fc)] + - ^ Ini?-/ (g) 



(106) 



The first two terms in Eq. ()106|) need to be expanded to the same way as in Eq. ((69|) . Keeping only the term, 
order A^. In the first, logarithmic term, we proceed in we obtain 



T^ln(-G„ 

k,(T 



{k)/W) = 



duj„i 



vpd (k — hp) 
A [vp {k - kp) - iiOmY 



duJr, 



A 



27r + A2 ' 



(107) 



r 



where Cjm = ojm — (wm) ■ The integral is controlled 
by frequencies O (A) , where the self-energy is small. Ne- 
glecting the self-energy, we arrive at the free-fermion-like 
result 



Tj2\n{-G-'ik)/W)^-- 



(108) 



k,(T 



Expanding the second, S^G term in Eq. (fT06l) to order 
A^ and integrating over k — kp, we obtain 



^ Kn) , ,^,2 - (109) 



2n 



A2 



from from energies of order W ^ where — S^(fc). How- 
ever, the low-energy E^(tj) is present in the last, hiD^^ 
term in (jllOp . which gives a nonanalytic contribution to 
to S. 

Minimizing Eq. (|110p with respect to A, and differenti- 
ating with respect to , we obtain the spin susceptibility 
as a function of T and H . Note in passing that the 
approach based on differentiation of the thermodynamic 
potential is completely equivalent to the diagrammatic 
evaluation of the linear susceptibility x(T, 7J = 0), used 
in earlier work,^'^°'^^ and also generates diagrams for 
the nonlinear susceptibility xiT^H). We illustrate this 
point in Appendix IB] 



The frequency integral is of order E^ {'-Om ^ A) / A, which 
vanishes in the limit A oo. 

We see, therefore, that S reduces to the sum of the 
free-fermion-like contribution (up to a renormalization 
of the Zeeman energy), the Hartree interaction, and the 
lnD|| term from the boson part Sb, 



1 



MA 



M2 



.9s, - 



vA^ 

~ ' 2 ' Av 



(110) 



D. 



Nonanalytic terms in the thermodynamic 
potential 



In this Section, we use Eq. (|110p to derive the 
nonanalytic terms in the thermodynamic potential in 
the vicinity of the critical point. We will see how the 
nonanalytic terms change in the non-FL regime. 



Comparing this expression with (|75)l and ([75]) . we see 
that the fermionic self-energy S^(cj) does not affect A^, 
MA and terms in S. This agrees with our earlier re- 
sult that PL renormalizations of these three terms come 



1. Away from criticality 

First, we discuss the FL regime, where E^ — ict;,„A. 
The logarithmic term in Eg. ljllOp reads 
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where = r2m(l + A) and 6 = 1 + > 0. 

Expanding the integrand of Eq. in 1 /q and evalu- 
ating the g— integral to logarithmic accuracy, we find that 
the {aq)'^ term under the logarithm can be neglected, so 
that Eq. reduces to 



"log — 



1 



(112) 

where ilg = <^ (j), — iA^ . Converting the Matsubara 
sum into a contour integral, we obtain 

where and are retarded functions of frequency, ob- 
tained via analytic continuation of the Matsubara func- 



tions 



/2 = [n-iA] In 



S\ n-iA 



^.2 



(.113) 



Performing the analytic continuation, we find 



;n (n + a) 



Im/2« = 77 f^ + A sgn 



Assembling the two parts, we obtain 



S + A] + g^n 



(114) 



"log 



dn coth ( — 



rn (^n + a) 



sgn 



6{n + A]+g''n 



(115) 



r 



In what follows, we choose A > without a loss of gen- 
erality. Because of the sign-functions, the integral in 
Eq. I|115p is confined to the interval 1/Ai > > A, 
where A 1 = l+g'^/{Xd). Since A oc l/\/J, we have Ai » 1 
close enough to the QCP. 

The thermodynamic potential Siog(T, A) at finite A 
and T is then given by 



5iog (t,a) 



2XT 



(116) 

This equation parameterizes 5iog(r, A) as a scaling func- 
tion of A/2XT. At low temperatures, T <C A/2A, one can 



replace coth(a;) by 1 -I- 2 exp(— 2a;) and extends the lower 
limit of the integral to zero. This gives 



Slog (t « A/A 



247ri}|,A 



O cxp - 




(117) 



In the opposite limit of high temperatures, T ^ A/2 A, 
one expands cotha; in series as cothx = -\- x/3 — 
x^/45 + .... The integrand is now logarithmically di- 
vergent, and the lower limit of the integral is relevant. 
Performing elementary integration, we obtain 



"log 



r> A/A =- 



InAi - . 



47rw|, 



A^T 



A4 



5767r vlX^T 



A6 



1728007r 4A4T3 



(118) 



Next, we obtain the thermodynamic potential ^ as a of the relation d^{M, A,T)/dA = 0, we obtain 
function of magnetization. Substituting (|117[) and (|118l) 
into (jllOp and expressing A in terms of M with the help 
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M2 

1^7 



(l+5s,o) 



|M| = 



1 



M2 

:(r»M/z.A) = — 



(l+5.,o)- 



T|ln(l 



(119a) 
(119b) 



The term in Eq. (jll9b|) determines the temperature 
dependence of the susceptibiUty 



1 /92^ 



/i| 2^1 



r|ln(l+g,,o)| 



(120) 

We can now compare Eqs. (|119ap and (|119b[) with Eqs. 
((79)) and ([SO]) , keeping in mind that one should set Wi? = 
Wi? in regular terms in Eqs. ()79I80|) as we measure the 
velocity renormalization with respect to its value at the 
upper cutoff for the spin-fermion model. We see that 
Eq. (|119a|) differs from Eq. (l79l) b y a factor of 1/A 
in the \M\^ term, while Eq. (|119bp differs from Eq. 
([50)1 by a factor of l/A^ in the M"* term. The factors 
are precisely vf/vf and {vf/vfY, respectively, in the 
spin-fermion model. This confirms our assertion that the 
nonanalytic term and its finite T equivalent M'^ /T 
are renormalized by fermions with the energies of order 
A = M/v. 

Equations (|119alll9b[) are valid in the PL regime, 
where max{T, M/i^A} ^ ojq/X^. At higher energies, the 
non-FL renormalization of the effective mass affects the 
functional form of the nonanalytic terms. This regime is 
considered in the next Section. 



2. At criticality 



propagator of spin fluctuations [see Eq. (jlOlj . Also, for 
energies in between ujo/^^ and wq, one can neglect the 
bare boson frequency compared to the self-energy, so that 
the logarithmic term in the expression for the thermody- 
namic potential becomes 



dqq 
2tt 



In 



v 



(121 



where Q — —icE^{flm,T), c w 1.2 (see [see Eq. l92b| . and 
the fermion self-energy is the sum of the quantum and 
classical parts, given by Eqs. ([55)1 and ([55)1 . respectively. 



a. T = 0. We consider first the T — case, 
when the Matsubara sum can be replaced by an in- 
tegral and the self-energy is purely quantum: = 



CLO, 



1/3 



sign(f]„)|rir 



12/3 



Introducing new variables 



- \ 2 

VFq/A] 



|A|3/2 



(122) 



The main difference between the non-FL- and FL 
regimes is the form of the self-energy, which enters the 



we rewrite Eq. (|12ip as 



c3/2|A|7/2 



dx I dyKe In 



VF \uja 



1/2 



(123) 



For a sufficiently small A, the first term under the log- 
arithm in Eq. (|123p can be neglected compared to the 
second one. We then obtain the field-dependent part of 
i^iog as 



. ,1/2 



(124) 



where z is the universal, i.e., cutoff- independent, part of 
the integral 



Re 



dx 



dy In [y 



2/3 _ 



if +X 



(125) 



This integral can be evaluated exactly and its universal 
part is equal to 



8\/27r 



(126) 
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so that 



^/2c3/2 |Ar/2 



357r 1/2- 



(127) 



Using Eq. ([57]) for cjq, expressing A via M, and adding 
the |Mp term, we obtain for the thermodynamic poten- 
tial as a function of magnetization M 



S(T = 0,M).— (1- 



.9.s.oj 



|M| 



7/2 



1/5/2^^/2 

where the energy scale is defined as^ 

4/3 



g 

afcc 



with 



-1 2/3 



hM^, (128) 



(129) 



5.22. 



(130) 



We see that S(r = 0, M) is still nonanalytic at QCP, 
but the leading nonanalyticity becomes |M|^/2 instead 
of |Mp in the FL-regime. Still, jMj''/^ is larger then 
the next-to-leading analytic term (M^). Correspond- 
ingly, the non-linear susceptibility scales with the mag- 
netic field as X cx -t-jHp/^. This scaling is dual to the 
|g|3/2 foj-m of the susceptibility at finite 

The crossover between the FL and non-FL forms of 
the thermodynamic potential [Eqs. ()119a|l and ()128|) ] 
occurs at the same energy where crosses over between 
the FL and non-FL forms, i.e., S is given by Eq. (|128[) 
for |Af| > u^vj}? and by Eq. (|119ap otherwise. In 
both cases, the nonanalytic terms are negative, which 
implies that a ferromagnetic quantum-critical point is 
intrinsically unstable against a first- order transition. 

b. Finite T. The form of the thermodynamic po- 
tential at finite temperatures depends on whether the 
temperature is above or below the scale Tqc, separating 
the regimes where the contributions from either finite or 
zero boson Matsubara frequencies dominate (T > Tqc 
and T < Tqc, respectively). In both cases, the non- 
analytic, |M|^/2 term is replaced by a regular, M'^ one; 
however, the temperature dependence of the prefactor of 
the is different in the two regimes. To see this, one 
can neglect the {aqf' term in Eq. (|12ip . integrate over 
expand the resulting expression to order A*, and convert 
the Matsubara sum into a contour integral. Following 
the same steps that led us to Eq. (jllSp away from the 
QCP, we expand coth(rj/2T) as 2T/n + 0/6T and keep 
only the second term in this expansion, which determines 
the coefficient of the A* term 



"log 



A4 
T 



1]2 



(131) 



For T > Tqc, f2 oc fl^/s ^nd the integral in Eq. 
scales as T^/'^. Correspondingly, the term in S is 



-M4/T1/3. For T < Tqc, ^ cx (T/|lnT|)^/^ and the 
AT' term behaves as -M^j lnT|. 

In addition, the T-dependence of the term, which 
determines the T-dependence of the (inverse) spin sus- 
ceptibility, changes from — T| ln(l -1-58,0) I at small T to 
— T| lnT| at higher T. The crossover occurs at T Ti, 
where Ti cx e^/A^ (Ref. iSlI ). Since the fermion self- 
energy enters the A^ term only under the logarithm, the 
difference between the regimes T > Tqc and Ti < T < 
Tqc is only in the numerical prefactor of the T|lnT| 
term. As it was pointed out in Ref. [Sll . the negative 
T| lnT| dependence dominates over the T-dependence of 
X within the HMM theory, which is given by the square 
of the thermal correlation length x^mm ^ ^^^C^) 
T/l lnT| and is weaker by a factor of (InT)^. 



E. Phase diagram of a ferromagnetic quantum 
phase transition 

In this Section, we analyze two possible scenarios for a 
ferromagnetic quantum phase transition in 2D, namely, 
the breakdown of a continuous transition and the spiral 
instability of a uniform magnetic state. 



1. First-order phase transition 

First, we assume that finite-g fiuctuations of the or- 
der parameter are negligible, and analyze a potential in- 
stability of a continuous second-order phase transition. 
Usually, the effect of nonanalyticities in the free energy 
on phase transitions is described in terms of the Landau- 
Ginzburg functional of the type given by Eq. (ISip . where 
the prefactors of regular (quartic and higher order) terms 
are assumed to be determined by fermions with energies 
of order of the bandwidth.- Phenomenological nature of 
these prefactors makes it difficult to make specific pre- 
dictions in this approach. Here, we will follow a different 
approach and calculate the entire thermodynamic poten- 
tial in a model with a long-range exchange interaction 
of radius a ^ kp^. The downside of this approach is 
the choice of a particular model. The upside is that not 
only nonanalytic but also analytic (Af^jM'*, etc) terms 
can be found explicitly, and therefore-at least within this 
model-one can make certain predictions about the nature 
of the phase transition. 

We will show that, in this particular model, the char- 
acteristic energy scale corresponding to the first-order 
phase transition is larger than the scale for the spiral in- 
stability, and falls into the regime where the mass renor- 
malization is small by a factor l/{akF) ^ 1. For a case 
when the interaction is short-ranged {akp ^ 1), both 
the first order transition and the spiral instability occur 
in the strong coupling, critical regime, and which one oc- 
curs first depends on the (unknown) prefactor of a regular 
term. 
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For the case afc^? ;:!> 1, we first assume and then verify 
that the mass renormalization factor A in Eq. (|128p can 
be set to unity. We will also see that the jump of spin 
polarization at the first-order phase transition, while still 
small compared to its maximum value (unity), is large 



enough so that one should analyze the full thermody- 
namic potential S(M, T) rather than its expansion up to 
order M**. Keeping in mind these two points, we write 
the T = thermodynamic potential as a function of M 
as 



S(M) 



(2^ 



In {aq) 



iMjvf + [vpqf 



(132) 



In Eq. (|132p . we set the coupling constant g to unity 
and neglected (5 under the logarithm; we will see that 
typical values of (aq)^ are much larger than (5. Equation 
(|132p can be reduced to a dimensionless form by rescaling 
X — {aqf' , y — flmCi/vFj C = Ma/vvp, and E {(,) — 
4a^S (A/) /i^Vp. An expansion of the logarithmic part 
starts with the term of order i^M'^ /v^ akp- We absorb 
this term into a renormalization of 6 in SM"^ / [Av), so in 
all formulas beyond this point 5 is already a renormalized 
parameter. Subtracting off the M-independent part, we 
obtain for the dimensionless thermodynamic potential in 
these variables 



aki 



-v{C) 



(133) 



where 



v{C) = - 



and 



1 



/"OO 

dy / dx In 

1 Jo 



(134) 



d = - [ dyy 

Jo Jo 



dx^ 5 ^ « 0.63. 

.iy^ + xf{x,/^x + yy 



(135) 



For small C, the expansion of V (C) starts with a nonan- 
alytic term: V {( <^ I) = |CI /3 + O (C^). Substituting 
the first, I CI /3, term into (|133p . we reproduce the non- 
analytic, cubic term in Eq. (|119ap . At larger function 
V {() goes through a maximum and falls off at C ^ 1 (cf. 
Fig. [21). 

The first-order phase transition occurs when the mini- 
mum in the thermodynamic potential at finite magnetiza- 
tion approaches zero (cf. Fig. ([6]), i.e, when the following 
two conditions are satisfied simultaneously 



£;(Ccr) = Oand £;'(Ccr) = 0. 



This yields 



2(5. 



ViC„) /akp 
2 



akp 



CV (Ccr) 

akp 



0. 



(136) 

(137a) 
(137b) 



Substituting p37ap into (jl37bp . we see that the jump of 
magnetization corresponds to an extremum of V (C) 



V (Cc 



0. 



(138) 



In our case, this extremum is a maximum. The first-order 
phase transition occurs when the line dc^o-kp touches the 
maximum of V (C), as shown in Fig.[71 Numerical calcu- 
lation gives 



Ccr ~ 1-11 and Set ~ 0.2l/akp. 



(139) 



Since ^cr ^ 1, a critical point cannot be determined 
by expanding ^^(C) in C s-nd keeping only a few first 
terms. Coming back to dimensional variables, we see 
that the jump of spin polarization at the transition 
Mcr/n ^ 1/akp ^ 1. The effective Zeeman splitting 
at the transition is Acr ^ ep/akp <ti ep. 

We now verify whether neglecting of the mass renor- 
malization was permissible. The critical distance to the 
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FIG. 6: (color online) Schematic behavior of the thermody- 
namic potential H as a function of magnetization M. Upper 
curve: the first-order phase transition has not been reached 
yet but the system exhibits a metamagnetic transition in finite 
field. Lower curve: the first-order phase transition occurs. 



FIG. 8: (color online) Phase diagram of a ferromagnetic phase 
transition. The lines of second- and first-order transitions 
(solid and dashed lines, respectively) are separated by a tri- 
critical point. 





FIG. 7: (color online) Function V{C,) defined by Eq. (IT34|l . 
The first-order phase transition occurs when the maximum of 
V reaches the value of 5 x [akp). 




QCP 5ci ^ l/akp corresponds to A ^ l/^/akp <C 1. 
Therefore, mass renormalization is, indeed, irrelevant. 
Another way to see this is to notice that the Zeeman split- 
ting at the transition Acr ~ ep/akp is parametricaUy 
larger than the characteristic energy ^ / {o-kp^ , 
separating the regimes of weak and strong mass renor- 
malization. 

At finite temperatures, the first-order transition occurs 
along the dashed line in Fig. [51 which starts at (5 = 5rr at 
terminates at a tricritical point {5tc,T = Ttc) (Ref. |62|. 
[63l ). For T > Ttc, the transition is second-order (solid hne 
in Fig. [5]). The tricritical temperature can be estimated 
from the condition that the M'^ / (v'^VpT) term, which 
replaces the |Mp one at finite T, becomes comparable to 
the regular bM'^ term. In our model, b ~ a/v'^Vp. This 
gives Ttc ^ vp/a ^ ^cr- To find the numerical prefactor, 
one needs to evaluate the entire thermodynamic potential 
at finite T but we are not going to dwell on it here. 



2. Spiral magnetic phase 

We now turn to another scenario of a phase tran- 
sition, which is possible due to a nonanalytic be- 
havior of the spin susceptibility as a function of the 
wavenumberiiiSii^i^ Away from the QCP, xil) scales 
as \q\] near the QCP, the non-FL mass renormalization 
changes this scaling to \q\'^^^. As in the previous Section, 
it will turn out that for akp ^ 1, the instability occurs 
before the system enters into the l^l^^^ regime. There- 
fore, we start with a FL form of xio) derived in Ref. [sil . 
In the long-range interaction model. 



^-\q)^S+iaqy -—(l + X)^. 



(140) 



Stt ^ ' ' kp 

Equation (|140p is valid for a moderate mass renormal- 
ization: (1 + A) ^ 1/VS- A spiral instability as g = <7cr 
occurs when conditions = ^ind dx~^/d\q\ = are 
satisfied simultaneously. This gives 

2 1 + A 

9cr = — 1— T (141) 



4 (1 + A)2 
97r2 kW ■ 



(142) 



Since A depends on S itself, as specified by Eq. (I57)) . 
Eq. (|142p represents an equation for 61^.. Solving this 
equation, we find that the spiral instability occurs at 



Sl= I ^ + J^^ + ^ I ^^»0.52/(afcF) . 

(143) 

At this (5, A - 1 and (1 + A)^^ - J < 1, so that Eq. (fMO]) 
is, indeed, applicable. 

Comparing the critical values of S for the spiral insta- 
bility and the first-order phase transition [cf. Eq. (jl39p ]. 
we see that 




(Sny 27r / (akpY 



(144) 
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FIG. 9: (color online) Re-entrant phase diagram of a second- 
order ferromagnetic phase transition. The two curves corre- 
spond to two solutions of Eq. (|145p . 



FIG. 10: Renormalization of vertex V{k,p\p,k) (hatched) in 
the Cooper channel. Dashed lines are irreducible Cooper am- 
plitudes. 




This means that the first-order transition occurs before 
the spiral instability. 



3. Re-entrant second-order phase transition 

Another interesting consequence of the nonanalytic be- 
havior of the susceptibility is a re-entrant second-order 
phase transition. This effect occurs due to a nonanalytic 
temperature dependence of x- Suppose that we are still 
far away from the tricritical point, so that the nonanalyt- 
icity of S as a function of M does not play a role. Taking 
into account the regular T dependence of x, which arises 
from energies of order of the bandwidth, we have from 

(inni) 

X-\T) = xo [5 + [T/nf - [T/ep) H\5\-')\ , (145) 

where To ^ W . The critical temperature of the 
second-order transition is determined from the condition 
X~^{Tc) = 0. In the absence of the nonanalytic term, 
the solution exists only for 5 < dX Tc — Tq\/^ and the 
transition line Tc{S) has a negative curvature for negative 
5. Due to the nonanalytic term, there exist two solutions 
of x^H^c) = for (5 > (see Fig. ^. One of them, 
Tci, vanishes as S/ln[d^^) at (5 — > 0, while the second 
behaves as Tc2 ^ In ((5^^). The two branches match at 
some (positive) value of ^ = S^. As a result, the phase 
transition occurs at 5 > 0, and the phase diagram ex- 
hibits a re-entrant behavior. Depending on whether the 
tricritical point is to the right or to the left of the reen- 
trant point, the phase diagram has the forms shown in 

Fig. m 
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F. Cooper channel near the quantum critical point 

We now address the issue of Cooper renormalization 
near a ferromagnetic QCP. We recall that away from the 
QCP, the backscattering amplitude vanishes at T as 
1/ ln[max{T, A}] due to singular renormalizations in the 
Cooper channel. Consequently, the backscattering con- 
tribution to the spin susceptibility is reduced by a factor 
of 1/ ln^[max{r. A}] as compared to non-backscattering 
contributions. The question is to what extent the Cooper 
renormalization affects the susceptibility near the QCP. 
This is an important issue because the sign of x(2^: 
determines whether the transition becomes first order or 
remain continuous. In the preceding discussion of the 
phase transition, we approximated the scattering ampli- 
tude by a single component fsfl, which diverges at the 
QCP, and neglected all other components. This is cer- 
tainly inconsistent with the vanishing of /s(7r) at T = 0, 
as it is the sum of all components that must vanish. We 
now include Cooper renormalization into consideration. 
For definiteness, we consider Sx{T,H = 0). 

First, we show that the logarithmical renormalization 
of /s(7r) starts below some temperature which becomes 
progressively smaller as the system approaches the QCP 
and \fs.o\ increases. To estimate this scale, we consider 
a vertex T{k,p;p, k). When k and p are projected onto 
the Fermi surface and corresponding frequencies are set 
to zero, r coincides with the scattering amplitude fs{0), 
which depends on the angle 9 between k and p. The 
backscattering amplitude corresponds to = tt. This 
vertex can be expressed as a series of diagrams in the 
Cooper channel (cf. FigfTU)). The dashed lines are irre- 
ducible amplitudes 7c (^) which, in general, depend on 
the scattering angle, e.g., on the angle 6' between k and 
1 in second term of the series. Integrating over the mag- 
nitude of 1 and summing over the internal fermion fre- 
quencies, we obtain for the Cooper bubble 



n,(0)(xiln^-, 
A cosO/2 



(146) 



where, we remind, A ^ 1 is the mass renormalization 
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factor. As a result, the series for fs{0) reads 



nent fs^i modelled as 



1 



A cos(e'/2) A 



cos(6i/2) 
(147) 

where 7^'' are the angular convolutions of two and 
more irreducible amplitudes; for example, {0) = 
J d4>^c{.4')lc{Q — etc. Since irreducible amplitudes 
7c(0) and ^^^\9) are regular functions of 9 near tt, they 
can be replaced by constants 7c (tt) and 7^^' (tt) for 9 x tt. 
Near 9 = tt, fs{9) then becomes 



fsiO) 



, £2 1 

A cos(0/2) 



C3 , 2 



cos(6i/2) 

= fsfi + Y.f^^r.cos{n9). (148) 



n>0 



The Cooper bubbles in Eq. (|148p have regular expan- 
sions over angular harmonics; for example In (cos 0/2) = 
— In2-|-^^^p(— )"'+^ cos(n0)/2n, etc. Equating the pref- 
actors of the cos(n0) terms in the first and second lines 
of Eq. (fT48l) . we obtain 



/s,o = ci + y In 2 



A 2n 



(149) 



According to our assumption of a Pomeranchuk-like in- 
stability, |/s,o| > 1 while |/s,„>o| ^ 1- Therefore, 
ci = fs,o + 0(1). Dropping terms of order unity, the 
first line of Eq. (|148p can be rewritten as 



fs{9) = fsfi + ^ In + ^ In^ 

A cos(6^/2) V cos(y/2) 



(150) 

The angular-dependent terms in Eq. (|150p are of order 
unity for ~ 1. As 9 approaches tt, they increase and 
cancel the large, angular- independent term, fs^. This is 
how the vanishing of /s(7r) happens. However, the can- 
cellation occurs only for angles exponentially close to tt: 
TT — ~ exp {—\fsfi\X/ci) . Associating typical it — (p with 
typical T/ep, we conclude that the characteristic tem- 
perature, below which Coooper renormalization becomes 
crucial, is exponentially small in |/s,o|A ~ (1 + gs,o)^^^^- 

These exponentially small temperatures are of no rel- 
evance for the preceding discussion, as they are lower 
that the temperature of the tricrical point, at which the 
second-order phase transition turns into a first-order one. 

Still, the very fact that Cooper renormalizations bring 
about higher harmonics of fs{(t>), which scale as Inep/T 
requires some attention because the divergent harmonic 
of the scattering amplitude , fs,o, enters 5x{T) only log- 
arithmically [see Eq. (pS] , and the presence of extra In T 
terms may affect the temperature dependence of x- 

To address this issue, we use a simplified model, in 
which partial components of fs{9) with n > are ab- 
sorbed into an effective temperature-dependent compo- 



/s,l — 



In 



T 



1 



A|/,,ol 



In^ 



(151) 



where t ^ 1 is a constant. This interpolation for- 
mula reproduces correctly the limiting forms of fs,i- 
Indeed, if In (ep/T) /A|/s,o| is small, fs,i scales as 
(l/A)lneF/T, consistent with the perturbation theory. 
If ln(eF/r)/A|/s^o| is large, fs,i approaches /s,o, consis- 
tent with the vanishing of /s(7r) = fsfi — fs,i sA. T — 0. 

We now use Eqs. ((50)) and (|5ip for the susceptibility 
and compare the contributions from g and f^ i. Before 
substituting fs^i from Eq. (|15ip into Eq. ([50)) for the 
susceptibility, we need to account for the mass renormal- 
ization in the expressions for the functions ^i(/s,o) and 
-F2(/s,o), which amounts to replacing fs^ by /s,o/A in the 
arguments of Fi and F2. Since i^i(/s,o) scales as In |/s,o|j 
it is not affected by mass renormalization. On the other 
hand, the asymptotic behavior of F2{fs,o) ~ 2/s^o in 
the absence of mass renormalization must be replaced 
by2/,,o/A. 

Using the modified expressions for Fi and F2 in 
Eq. (|50)) and the perturbative form /c,i ~ (c/A) Inep/T, 
valid for all but exponentially small T, we find 



Sx{T,H = 0) = x'o 



2D 



Ms 



T 



(152) 

This is the key result. We see that Cooper renormaliza- 
tion generates a InT correction to our earlier result for 
x(r) in the isotropic scattering model. However, at all 
but exponentially small T the correction is small in 1/ A^ 
and can be safely neglected. We conclude, therefore, that 
the sign and magnitude of x(T, H = 0) near the QCP is 
not affected by renormalizations in the Cooper channel. 

The same consideration and conclusion hold also for 
X{H.T = 0). 



IV. MAGNETIC RESPONSE IN 3D 

In this Section, we discuss the nonanalytic behav- 
ior of the spin susceptibility in 3D. The behavior of 
SxiT\qJI = 0) has been considered in Refs. HIIllIll, 
[T9II34I35II63I65I . As the Kohn anomaly is logarithmic 
in 3D, the corresponding nonanalyticities are also only 
logarithmic: Sx{T,q,H = 0) cx q^lnS, where £ — 
meix{vF\q\,T}, and Sx(T = 0,q = 0,H) (x H^ln\H\. 
It was argued in Refs. IjIsI that these non-analyticities 
give rise to tricritical points, observed in 3D itinerant 
ferromagnets.^ However, the relative weakness of nonan- 
alytic terms in 3D raises a concern whether the tricritical 
behavior is indeed caused by the many-body physics or 
is due to some single-particle effects, such as features in 
the density of states^ or magnetoelasticity^i. In this 
Section, we show that logarithmic nonanalyticities are 
weakened even further by the mass renormalization in the 
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vicinity of a 3D FM QCP, to q^ln\ln£\ and In | In | F| | 
forms, respectively . 

Another peculiarity of the 3D case is that the T- 
dependence of }(iT, H = 0,q ~ 0) is analytic: Sx{T) oc 
C3dT^ (Refs. 34.35. ) . For a Galilean-invariant system 
(with a quadratic fermion dispersion), such as He'^, the 
coefficient c^d was shown to be negative in Ref . [3^. i.e., 
the T-dependence of x has the same sign as in a Fermi 
gas. We reanalyze this result in Sec. II V 2l and show that 
the magnitude and sign of c^u are nonuniversal, i.e., they 
depends on details of the fermion dispersion. This may 
explain while c^d is positive in He'^ fRef. [69h but negative 
in some exchange-enhanced paramagnetsi^ 



1. Magnetic-field dependence of the susceptibility 

We first discuss the nonanalytic behavior of the 
susceptibility at T = in the perturbation theory. 

a. Perturbation theory To obtain the \n\H\ field 
dependence of the spin susceptibility in 3D at T = 0, it 
suffices to expand the integral for S to order and ver- 
ify that the prefactor diverges logarithmically. Cutting 
off the singularity at the scale set by H, one arrives at 
the In \H\ behavior of S and, therefore, the In \H\ 
behavior of Sx- 

At second order, S is still given by Eq. (O, where H^i 
is now the 3D polarization bubble. For simplicity, we 
assume that the interaction is local, i.e., U{q) ~ const = 
U. 

Similarly to the 2D case, IIi^ can be separated into the 
static and dynamic parts as 

-I- vpq + A 

Hn i^m, q) = In — — - 

(153) 

Here and till the end of this Section ly — rnkp/^TT^ is 
density of states at the Fermi energy in 3D. Expanding 
the square of H-^ to order A**, we obtain 



terms, and 



)(2) 



(x) = X- 



2 arctan(a;~i)x [x^ - l] + x"^ + 4/3 



1) 



(155) 

Substituting Eq. (|154p into Eq. ([7]) and integrating over 
firm we find that the g-integral is indeed logarithmic. 
Cutting off the logarithmic singularity by the field, we 
obtain 



S2(T = 0,ff) 



dq 



{2tt)^vf^ J\a\/ 
A^ 



487r2 vl 



In 



A 



dy Q (y) 



(156) 



where u = Uu and we used that J^^dxQ'^^^ (x) = tt/6. 
Consequently, 



Sx 



(2) 



44 A 



(157) 



where Xo^ = rnkpfi^/TT^ is the spin susceptibility of a 
3D Fermi gas. We see that that the sign of the field- 
dependence of the second-order contribution is metam- 
agnetic, i.e., y increases with H. This conclusion contra- 
dicts to Ref. [to, where x was found to decrease with H as 
(without a logarithmic factor) . Notice also that typi- 
cal r2„j and q in Eq. (|156p are such that fl„i ^ vpQ S> A. 
Therefore, even the second-order contribution in 3D does 
not arise entirely from backscattering processes. 

Higher-order terms modify the second-order result by 
renormalizing static vertices in the skeleton diagram with 
two dynamical bubbles, and add new processes involving 
larger number of dynamic bubbles. For example, the 
nonanalytic part of diagram e in Fig. [2] contains the cube 
of the dynamic bubble. Expanding H^^ to order A"*, we 
obtain 



H^, (f]„, g) = . . . + ^2 fA.) Q(2) + . . . (154) 

where dots stand for A-independent and regular A^ where 



H^^ inra,q) = 



VFq 



\VFq 



(158) 



, , ,3x^ (x^ -I- 1) arctan^(a; M + x (4 + 3a;^) arctan(a; - 3a;^ - 1 
J*-"*^ (x) = x^ ^ 



(a;2 + 1)4 



r 



(159) 



Substituting this expansion into Eq. (j33[) , and performing Eq. (jl56[) , we obtain 
the integrations over Sim and q in the same way as in 



8x 



(3e) 



44 A 



(160) 
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where we used that J^^dxZ {x) = 7r(7r^ — 8)/16. As 
it also the case in 2D, the signs of the second- and 
third-order contributions are opposite. 

b. Quantum- critical behavior in 3D The vicinity of 
a ferromagnetic QCP in 3D can be analyzed within the 
Eliashberg theory, in the same way it was done in Sec. lIIII 



for the 2D case. Since our goal here is to obtain only the 
qualitative behavior of X; we will not repeat the manip- 
ulations within the spin-fermion model but just consider 
the RPA for S, which reproduces correctly both the weak- 
coupling and quantum-critical limits. Summing up the 
RPA-series and replacing u by —gs,o we obtain 



S(T = 0,i7) = ^ 



dflr, 



J 



q^dq In 



.5 + 



irir 



In 



vpq 



2vFq i^r. 



vpq 



A 



(161) 



where F = —gs.o/{l + 5s. o) aud dots stand for regular 
terms. Expanding the integrand to fourth order in A, 
and integrating over Q^n and we obtain 



where 



■(r = o,i7)^-^|-z(r)in 



A 



Z{T) 



w/2 



(162) 



cos 2(j) 11 cos^ 0—2 cos^ 



652 



53 



4S'4 



(163) 



and S ^ T-^ - 0cot0. For small F, Z{T) « 
and Eq. pM]) reduces to Eq. (fT56)) . At the QCP, F ^ 
00, and K{r) « 0.032. Substituting this into (flB^ and 
differentiating twice over H, we obtain at the QCP 



This is a very weak nonanalytic dependence. 



(164) 



Equation (|164p is, however, incomplete as the RPA ne- 
glects mass renormalization which becomes singular near 
the QCP. Near a 3D QCP, the fermion self-energy be- 
haves as T,{ujjn) oc w„i|lnWm|, i.e., at small frequencies 
it is parametrically larger than a bare term in the 
fermion propagator. Reevaluating the polarization bub- 
ble for dressed fermions, we obtain the same expression 
as in Eq. (jl53p . except for under the logarithm is 
replaced by In \ ilm\- Substituting this into Eq. (|16ip . 
setting F = cxD, expanding again to order A"* and inte- 
grating first over n^n aud then over q, we obtain 



S2(r = O.i?) - A^ / - A'^lnllnlAI 



Consequently, 



A 



(165) 
(166) 



2. Temperature dependence of the spin susceptibility 



The temperature dependence of the spin susceptibil- 
ity in 3D, Sx{T) oc T^, was found perturbatively in 
Refs. [7II35I and in the paramagnon model by Beal-Monod 
et al. (ReflU). Here, we consider the temperature depen- 
dence of the spin susceptibility of a 3D FL near a ferro- 
magnetic QCP. As in Ref. Is^, we consider fermions with 
local interaction U near a ferromagnetic QCP, so that 
Uv = —gs,o ~ 1- However, in contrast with Ref. [3^ . we 
assume an isotropic but otherwise arbitrary energy spec- 
trum £fc. Correspondingly, the density of states, v (efe), is 
an arbitrary function of the energy. As it will be shown 
in this Section, this generalization leads to a possibility 
of reversing the sign of the T-dependence of x- 

The thermodynamic potential is given by Eq. (j65p : 
however, in contrast to the 2D case, the T-dependent 
part of the susceptibility comes not only from the ladder 
but also from the ring diagrams for S. Consequently, the 
(5S(T) term in Eq. (|65p contains both types of diagrams 
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FIG. 11: Diagrams for the thermodynamic potential that give 
rise to the temperature dependence of the spin susceptibihty 
in 3D. Top row: ladder diagrams; bottom row: ring diagrams. 








(cf. Fig.[TTl): SE = SEl + SEr, where 

5El = T^[ln(l + C/nu)] (167a) 

SEn - ir^[in(i-c/2ntni) + c/2nTnJi67b) 



and 



Il^,l{n„,, q) ^Tj^G^lik + 'z)GT.i(fc) 



(168) 



is the polarization bubble composed of fermions with the 
same spin. The second-order diagram can be considered 
as a part of either the ladder or the ring series; we asso- 
ciate it with the ladder series and subtract the second- 
order term from the ring sequence. 

The details of the calculation are presented in Ap- 
pendix [C] Here, we present the results for two different 
regimes near a QCP. The first one is the FL regime of 

smah T and finite l + gsM- (1 + gs,of^^ » T/ep ■ In this 
regime, SxiT) scales as T^, and the prefactor is the sum 
of the universal (u) and nonuniversal (n) contributions: 



T 

l + {Au+An) ( — 



(169) 



where x(0) — (1 + gs.o). The universal contribu- 

tion has the same origin as in 2D: it arises due to a long- 
range dynamic interaction between fermions mediated by 
particle-hole pairs. The prefactor of this contribution, 
Au, is given by 



1 



3 {i + gs,oy 



(170) 



The nonuniversal part comes from the regular terms in 
the expansion of the polarization bubble in the momen- 
tum and frequency. The magnitude and the sign of the 
nonuniversal part depend on details of the fermion dis- 
persion: 



Ar, 



7r2 (v'y - (5/3)z/i/" 



(l+5s,o) 



(171) 



where Vk = dsk/dk, rrif,^ = cPsk/dk^, v' = 
dv{£k)/d£k\ek=eF, = d^^i^k) / del\e^=ep ■ The coeffi- 
cient a determines the q dependence of the static bubble 
for free fermions in 3D 



n(an = 0,9) = -V 



12 Kkp 



(172) 



and is given by 



kl d 



d£k 







1 " 









2 d 

3 dsk 



(173) 

For a power-law dispersion, Sk = Ak'', we have a = 
(7 -I- l)/3. For a quadratic dispersion, 7 = 2, a = 1 and 



Ar^ 



1 



6 (l+5s.o)2 



(174) 



Comparing (|170|) and (|174[) . we see that An is positive 
and larger in magnitude than A^: A^jA^ = k, 
—4.93. This is the result found in Ref. '34. In this case, 
X decreases with T: 

X(T) = 



X(0) 



1 



l + 9s.O \<^F J 



(175) 



This case is perhaps most relevant for He'^. However, for 
more complex dispersions (relevant for anisotropic FLs in 
metals), the universal term can win over the nonuniversal 
one. Indeed, because a must be positive regardless of the 
dispersion (otherwise, a system of free fermions on a lat- 
tice would have a magnetic instability at finite (7), the sign 
of An is determined by the the sign of the combination 
{v'f - (5/3)z/i/", entering Eq. PTT|) . We see that A„ > 
as long as v" < 0. On the other hand, if the Fermi en- 
ergy is near the minimum in the density of states, where 
v' ~ {) and v" > 0^ then An < and x increases with T. 
Therefore, both types of the T-dependences are possible 
in metals. An increase of x(T) with T was observed in a 
number of strongly paramagnetic metals"^^. 

Next, we consider the quantum critical regime: 
T/ep ^ (l + gsfi)"^^^. Here both universal and non- 
universal terms scale as T^^^ (Ref. [68h . A detailed cal- 
culation shows that 



X-\T) = (2Ai» ' b [Bu + Bn) {T/ep) 



4/3 



(176) 



where B^ = — S/tt^ is the universal, low-energy contri- 
bution, Bn is the nonuniversal contribution whose form 
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depends on the details of the fermion dispersion, and 
b = 22/3(67r)i/3r(4/3)C(4/3)/35/2 « 0.871, with T{x) 
and (^(x) being the F- and C-functions, respectively. For 
the quadratic dispersion, i?„ = 1 > S/tt^ and the suscep- 
tibility is positive, i.e., the system is stable. It is possible, 
however, that B is smaller than Stt^ for a more complex 
dispersion, in which case the susceptibility is negative at 
the QCP. This implies that the transition line Tc{gsfl) 
may have "wrong" sign of the slope at small T, as in 
Fig. [HI and the system exhibits a reentrant ferromagnetic 
transition before reaching the QCP at T = 0. 

There is also a more fundamental reason for the re- 
entrant behavior. In the calculation that led us to 
Eq. (|176|) . we neglected the fermion self-energy, which 
scales as LOm In LUm at the QCP. Including the self-energy, 
we find that the universal, negative contribution to the 
inverse susceptibility [the — S/tt^ term in Eq. (|176p ] ac- 
quires an extra factor of (InT)^, while the non-universal 
contribution remains the same. As a result, the universal 
term wins at low enough T, and the system displays a 
re-entrant behavior even for the quadratic dispersion. 



V. CONCLUSIONS 

In this paper, we have considered nonanalytic spin re- 
sponse of an interacting Fermi system, both away from 
and near a ferromagnetic quantum critical point. Our 
two primary goals were i) to establish the sign of the 
magnetic-field and temperature dependences of the spin 
susceptibility in the Fermi-liquid regime and ii) ana- 
lyze the stability of the continuous ferromagnetic quan- 
tum critical point in 2D. We found that higher-order 
processes, involving more than one particle-hole pair, 
may reverse the anomalous (positive) sign of the single 
particle-hole pair contribution in a 2D Fermi liquid. A 
controllable calculation within a large-iV model shows 
that this effect is more important than Cooper renor- 
malizations of the backscattering amplitude considered 
in Refs. [TgHTsI . For a 3D Fermi liquid, we showed that 
the sign of the T-dependence of x is determined by the 
balance between universal and non-universal contribu- 
tions; while the former depends only on the density of 
states and Fermi velocity near the Fermi level, the latter 
is sensitive to the actual form of the fermion spectrum 
away from the Fermi energy. We found that whereas the 
sign of T and H dependences of x in the Fermi liquid 
regime depends on the detailed form of the interaction, 
the anomalous (positive) sign of these dependences is re- 
stored near a ferromagnetic quantum critical point. At 
the same time, the role of Cooper rcnormalizations is 
diminished even further near criticality. Analyzing dif- 
ferent mechanisms for a breakdown of second-order phase 
transition, we showed, within the model of a long-range 
exchange interaction, that the first-order phase transition 
preempts an instability towards a spiral magnetic phase. 
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APPENDIX A: SPIN SUSCEPTIBILITY OF A 2D 
COULOMB SYSTEM IN THE LARGE N LIMIT 



In this Appendix, we calculate the nonanalytic 
magnetic-field dependence of the spin susceptibility for 
a valley-degenerate 2D electron gas with Coulomb in- 
teraction. We assume that there are degenerate or- 
bital valleys, so that the total (spin x valley) degeneracy 
is = 2Ny ^ 1. For simplicity, we consider only the 
T = case. 

a. Ring diagrams Since each polarization bubble 
comes with a large factor of the leading contribution 
to the thermodynamic potential is given by a series of 
ring diagrams [diagrams a,b... in Fig. [2], which contain 
a maximum number of bubbles at each order. The sum 
of ring diagrams gives 



i-^uiq) (n^ 



Hi) 



(Al) 



where U{q) — 27Te'^/q is the bare Coulomb potential. The 
nonanalytic dependence of S is determined by q near 2kF- 
For q « 2kp, the second term under the logarithm in 
Eq. (|Aip is of order N{e'^ /kF)Tn ~ gN, where is the 
dimensionless coupling constant defined in Eq. (j6ip . For 
gN <C 1, one can expand the logarithm to second or- 
der in U{q). This reproduces a weak-coupling result of 
Ref. [111. In the opposite limit oi gN ^ 1 (strong cou- 
pling), the second term under the logarithm dominates, 
and the field-dependent part of ^n=oo reduces to 



iT^ln[l-(Pi.+PJ/2], 



(A2) 



where P|,| are the dynamic parts of the bubbles with the 
same spins, defined as 



nT.i(a, 



,q)^--il-P,,0- 



(A3) 



To evaluate Hjv^ooi we expand the logarithm in 
Eq. (jA2p in Taylor series 
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2 (^T + Pi 



-\ 3 



(A4) 



r 



This expansion can be viewed as a series of fictitious ring 
diagrams with each bubble being either P| or and the 
effective interaction being equal to unity. 
At r,ir = and g« 2fcF, 



where q — 
\q\vF, Eq. 




2fci? and |(7| 
reduces to 



< fcF (Ref. 0)- For 



P|_l instead of the full bubbles. Since P|,| do not contain 
constant terms, diagrams containing more than two P|,| 
do not give rise to a nonanalyticity. 

An explicit computation of S7v=oo from Eq. (jA7p is 
quite involved, because one needs to know the expressions 
for 2kp bubbles in a finite field. Fortunately, there is no 
need to perform such a computation because the nonan- 
alytic part of the second-order diagram for a short-range 
interaction 



^2 = -^c/^T^ n^n^ 



(A8) 



Pt - Pi 



{q/kp f^ ; for g > 0; 
^/2\n^\/2iF{kF\q\f^, forg< 0. 



One can easily verify that the field dependence of the 
^N=oo comes only from the products of with oppo- 
site spins. At second order, there is only one such a term: 
P|P|. The g- integral of P^Pi diverges logarithmically for 
q < 0. Finite magnetic field splits the Fermi momenta of 
spin-up and spin-down fcrmions. This, along with with 
finite f2m, regularizes the integral: 



- max{|i2„|,|A|} 



dq/\q\ 



= -02^1nmax{|a„|,|A|}. (A6) 

This is the same logarithmic behavior that we have 
already encountered in the perturbation theory [cf. 
Eq. (fTO|) ]-it leads to a | term in S. At third order, the 
nonanalytic dependence on A comes from PfPi- The q 
integral of P^P^ yields jfi^p / dg/|g3/2 ^ This 
leads to a |A|^/^ nonanalyticity in S, which is weaker 
than the | Ap one. All higher order contributions produce 
either analytic terms or nonanalytic terms with higher 
powers of A. Therefore, as long as only ring diagrams 
are considered, the only source of a |Ap nonanalyticity 
is the second order term in (IA4I). Hence, 



(A7) 



This situation is to be contrasted with the regular pertur- 
bation theory [cf. Sec. Ill A] , where ring diagrams are built 
from the full bubbles Hi.i (see IA3p . There, diagrams to 
all orders yield |Ap terms. For example, at third or- 
der one generates a |Ap term by keeping P|,| parts in 
two out of three bubbles and replacing the third one by 
its constant part (— m/27r). A strong-coupling expansion 
that we consider here involves ring diagrams built from 



contains the same combination of P|Pj. 
Therefore, 



as Eq. 



f . 



(A9) 



where u = vU . Using Eq. (jTU]) ^2 and differentiating 
twice with respect to the field, we obtain 



5XN=oo = ^ 



Xo ■ 



(AlO) 



Notice that this result could be obtained from the 
second-order ring diagram by replacing the wavy line by 
the screened Coulomb potential for ^ 00 



U 



2^e2 



27r 

771 A^' 



(All) 



The 2k p ring diagrams can be viewed as series of 
backscattering events, hence Eq. (jAlOp is the backscat- 
tering contribution to the susceptibility for a Coulomb 
interaction in the large- A^ limit. 

h. Cooper channel Out of the subleading, 1/A^ 
diagrams that modify the backscattering contribution 
to 5x^ the most important ones are the Cooper dia- 
grams which come with an additional logarithm, Lc — 
ln(Vl^/|A|). At low energies, when Lc > N, the Cooper 
diagrams become comparable to the ring ones. 

To third order, a Cooper diagram is diagram c in Fig. [2] 
with wavy line replaced by the screened Coulomb inter- 
action from Eq. (|A11[) . This diagram reads 



C \^ 

n) ■ 



(A12) 



where 11^'^ is a Cooper bubble composed of fermions with 
opposite spins 

n^l=Tj2G^ik + q)Gi{-k). (A13) 
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At T = 

n 



W 



27r 



m 
2^ 



f7™ + i A + y (a« + jA)2 + 

[Lc + Pfl] , (A14) 



where P^^ is a dimensionless function of two variables, 

2 



a; = rim/|A| and y 



pC 

Mi 



-Rein 



X + isgnA + y {x + isgnA)2 -)- y 



Substituting Eq. (IA14P into Eq. (|A12p . we obtain 
T 



(A15) 



12iV 



3Z/cP|j 



-3Lc [Pa%[Pi 



(A16) 



Cooper rcnormalization of the backscattering amplitude 
comes from the third term in Eq. (|A16p 



(A17) 



while other terms either give field-independent contribu- 
tions or generate 1 /N corrections to the prefactor of the 
|A| term. We now recall that the second-order diagram 
(diagram a in Fig. [2|) can be represented equivalently 
either in the particlc-holc or the particle-particle form. 
The latter reads 



Comparing Eqs. (IA17I1 and (|Al8j) . we find that 



As a result, 



Lc 
2N 



2L 



c . 



^2 \u=l 



N 



C,3 



2L^ 



■Sx 



(A18) 



(A19) 



(A20) 



where SxN=oa is given by Eq. (jAlOp . 

Higher order Cooper diagrams (diagram d in Fig. [2] 
and similar diagrams at higher orders) form a ladder. At 

any order n, we need to select only the L^~^ 



pc 

Mi 



term 



from the n"^ power of the Cooper bubble. The n*^ order 
Cooper diagram gives 



i-Y' 



N\ n{n - 1)- 




(A21) 



Re-expressing again TJ2„ 



pC 

Mi 



via :^2, we obtain for 



the n order Cooper contribution to the susceptibility 

Sxc,n = (-)" (n-l) ( j SxN=oo, for n > 3. 

(A22) 

Summing up all orders, we obtain 



^Xc = X! ^^'^''^ = 



n—3 



{l + Lc/Ny 



(A23) 

Combining this result with Eq. (lAlOp . we find for the 
total backscattering contribution Sxbs = Sxn=oc + Sxc 



xr 



Sep (1 + Lc/N) 



2 • 



(A24) 



c. 1/N corrections The 1/N corrections to 
Eq. (|A24p come from diagrams of third and higher or- 
ders in the screened interaction (jAlip . excluding the 
main logarithmic parts of Cooper diagrams. Some of 
these diagrams, e.g., diagrams b,d and / in Fig. [21 be- 
long to backscattering type and contribute either non- 
logarithmic corrections to the backscattering amplitude 
or subleading logarithmic corrections in the Cooper chan- 
nel. For example, fourth-order Cooper diagram d con- 
tains not only the Lq term, already taken into account 
when summing up the leading logarithmic terms, but also 
a subleading, Lc, term. These diagrams change the two 
terms in the denominator of (|A24p as 1 — > 1 -f 0(1/A^) 
and Lc/N Lc{l/N + 0{1/N^)). At low energies, i.e., 
for Lc ^ N, the entire backscattering contribution is 
then given by an asymptotic limit of Eq. (|A24p : 



i^Xbs = 



\A[xrN^ 
Sep 



7^2 



(A25) 



In addition, there are diagrams of different type, which 
do not participate in the rcnormalization of the backscat- 
tering amplitude. To leading order in this is dia- 
gram e in Fig. [2) This diagram has already been calcu- 
lated in Sec. Ill B 21 we just need to substitute u = 1/N 
into Eq. p6)) and multiply the result by (iV/2)^. Combin- 
ing this result with the backscattering contribution, we 
obtain the low-energy form of Sx in the large-A?^ model 
given by Eq. 



APPENDIX B: EQUIVALENCE BETWEEN THE 
THERMODYNAMIC AND KUBO-FORMULA 
APPROACHES FOR THE SUSCEPTIBILITY. 

In this work, we were using a thermodynamic ap- 
proach to spin susceptibility, in which one first evaluates 
the thermodynamic potential in finite magnetic field and 
then finds the susceptibility by differentiating the po- 
tential with respect to the field. Alternatively, one can 
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find the susceptibility via the Kubo formula for the spin- 
spin correlation function. Certainly, the results of the 
two approaches must be equivalent in the linear-response 
regime. However, when going beyond the lowest-order 
linear-response diagrams, it is quite easy to miss cer- 
tain diagrams which are generated automatically in the 
thermodynamic approach. Moreover, the diagrams for x 
generated by differentiating the thermodynamic poten- 
tial work even beyond the linear-response regime, i.e, for 
an arbitrary ratio of the temperature and the Zeeman 
energy (but for Zeeman energies still smaller than the 
Fermi one, when the magnetic-field dependence of the 
bare interaction can be neglected). 

We now show how the diagrams for x ^i'g generated 
within the RPA in the spin channel for a local bare in- 



teraction: U{q) = const = u/v. In this case, the field- 
dependent part of the thermodynamic potential is given 
by Eqs. (|65p and (|66|1 . The polarization operator is given 
by Eq. @ with the Green's functions from Eq. (|14p . 
Differentiating Eq. ([66)1 and using the identities 

^ = TMbG?, (B1) 

we obtain 

5x{T,H) = 5xi + 5x2. (B2) 

where 



5x1 = 2A|t2 ^ i^un () (^)^? (k + q)- ik)G^ (k + q)- Gt {k)Gl {k + q)] 
q k M\q) 

5X2 - Mlr^EEE n , f/n (aW [G'(k,)G,{k,+q) - G^{k,)Gl{k, + q)] 

q ki k2 ' U\'i)) 

X [G^(fc2)Gx(fc2 +q)- G^{k2)Gl{k2 + q)] . (B3) 



FIG. 12: Diagrams for the spin susceptibility generated by 
differentiating the thermodynamic potential with respect to 
the magnetic field. 




Algebraic expressions for 5xi and 5x2 are equivalent to 
diagrams in Fig. 1121 The first term in i^xi is equiva- 
lent to diagram a for the vertex correction insertions into 
the bare particle-hole bubble; the second and third terms 
are equivalent to diagrams and b and c, which are self- 
energy corrections. The contribution 5x2 is equivalent 
to " Aslamazov-Larkin" diagrams d and e, which contain 
two triads of fermion propagators, connected by the in- 
teraction lines. The wavy lines in these diagrams are spin 



propagators 1/(1 + UIii^i{q)), while the vertices contain 
Pauli matrices ctq/j. The same diagrams have been used 
in Ref. [Sll to analyze diagrammatically the momentum 
dependence of the spin susceptibility. 
APPENDIX C: TEMPERATURE DEPENDENCE 
OF THE SPIN SUSCEPTIBILITY IN 3D 

In this Appendix, we derive Eqs. (jlSQHlTSp . describing 
the temperature dependence of the spin susceptibility in 
a 3D FL near a ferromagnetic QCP. 

As we said in the main text, the field-dependent part 
of the thermodynamic potential is represented by a sum 
of the ladder and ring diagrams, which contain polariza- 
tion bubbles made of fermions both with the same and 
opposite spins: (nx.x)) and IIxx, respectively [cf. 

Eqs. (|167all67bp ]. In 2D, the nonanalytic, 0{T) behav- 
ior of the spin susceptibility is associated with the field 
dependence of the the dynamic part of IIxx • ^ 
dependence of x turns out to be analytic, and one needs 
to consider the field dependence of both static and dy- 
namic parts of IIxx, as well as the field dependence of IIx 
and Hx (cf. Eq. (fTBg]) ). 

We begin with IIxx- Expanding nxx(0,0) to order A^, 
we obtain 
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nu(o,o)= /-^/^GTK.fc)G,K.,fc) = i _ 

{2n) J ^TT A J-A/2<£fc-CF<A/2 (27r) 



24 



(CI) 



r 



We will also need a term in the expansion of a static ence SIl{q) — — Q,q) — Il{flm = 0,9^0) to order 
bubble in zero magnetic field. To obtain this term for an : 
arbitrary but isotropic dispersion, we expand the differ- 



J 



m(g) 



i2nf 



fa (£k+q) - / (£k) 

£k+q — Ek 



/o(ek) 

I 



(2^ 



(C2) 



where /o(£k) is the Fermi function and 

Ss = £k+q — Ek- 



For an isotropic system et = £fc and 



(C3) 



£k 



e'k qx + 



2fc 



(C4) 



r 



For the quadratic spectrum, the last term is equal to zero. n^: 
Integrating by parts, we find 



(Pk 



q^ d 



.2 J2 



-V — V — — + VOL 



fill I \2 " \ \(^l\'^ 



{2n) 



3 dei 



(C5) 



24 12fc|, 

ii^r?™ if7„j + vpq + A 

+ — In —. (C7j 

^vpq iVLm - Vpq + A 



Adding up these two results, we obtain 



n(0, g) ^-v + va 



Next, we consider and Hi, which depend on A only 
via the static parts 



i2ky 



(C6) 



where a is given by Eq. (|173p . 

Combining Eqs. (jCip and (|C6[) with the dynamic part 



Ht,; (0,0) = -i/ eF±A/2 



of IIii from Eq. ()153|) . we obtain a complete result for Using Eq. p67a|) for S^LiT, H), we obtain 
I 



(C8) 



. _i 1 d^SEL 



[fJiBl^Y 9A2 A=0 (^si^) 



in 



172 



/an 



Ti 



i + f/nti aA2 (i + f/ntx)' V aA 
I 



A=0 



(C9) 



The derivatives of IIii in Eq. (|C9p are found with the help of Eq. (jC7l 



in 



9A2 



2;/- 



9n 



n 



9A 



A=0 



A=0 



12 {ni^ + vW) 
n2 



^{ClOa) 



(ClOb) 
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The first term in Eq. ( IC9p gives a nonuniversal con- expanding H^i itself to order {q/kp)^ ■ Doing so, we find 
tribution, which can be obtained by expanding the pref- 
actor (l + C/IIii) ^ to first order in Q.m/vpq and also 



12 {^LBV 



-T 2^ 1 + gsfi 



2vFq 12fc|, 
I 



1 Ppil v" 



12x(0) a V (l + g«,o) 



(Cll) 



where x(0) = 2/x^i// (1 + .gs,o)- The Matsubara summa- der In |1 + .gs.ol/ (l+gs.o); which diverges weaker upon 



tion was performed using 



T^|f]™| = const - ^T^ 



approaching the QCP than the r^/(l + gs.o)^ contribu- 
tion in Eq. (jCll[) . Therefore, this contribution can be 
(C12) neglected. 



The universal part of the ladder contribution comes 
The second term in Eq. (jC9P gives a contribution of or- from the second term of Eq. (jC9p . We have 

I 



2 3 ■ 



(C13) 



The total ladder contribution is the sum of Eqs. (jClip The contribution to x from ring diagrams is obtained 

and (|C13p . in the same way, and the result is 



ni^Bvf ^1 1- 



i72ntnx (i-c/2ntnx) 
I 



A=0 



(C14) 



Using Eq. (jC8[) . we find that the second term in Eg. (|C14p where 11 is the polarization operator in the absence of 
vanishes, while the first one reduces to the field. Obviously, the entire ring contribution is non- 

universal. 



^« \{[iBV? V ^ ^ (2^)^ 



1 — f/2TT2 

^ J Near a ferromagnetic QCP, i.e, when g^.o ~ — 
(C15) Eq. (jUTSl) can be simplified to 



^Xr =-^, T7^ ^E l + .9^.o 



2vFq I2k'j, 



1 k^vl {v'f - vv" 



2 ,,„3 



(C16) 



Combining with Eq. (jClip for the nonuniversal part of the ladder contribution, we obtain Eq. (fT7T|) . 
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while the universal part of 5xl ^ gives Eq. I|170p . 
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